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I  INTRODUCTION 


This  report  describes  a  one-year  research  and  development  effort  to 
determine  the  feasibility  of  using  hand  geometry  and  other  features  of 
the  hand  as  bases  for  the  construction  of  an  automatic  personnel  authen¬ 
tication  system  for  access  control.  The  subsequent  sections  of  the 
report  describe  the  method  and  extent  of  data  collection,  the  method  of 
data  preparation  and  reduction,  the  data  analysis,  the  experimental 
authentication  results,  and  provide  brief  descriptions  of  the  computer 
programs  that  were  written  for  this  work. 

The  appendices  are  photocopies  of  listings  of  the  computer  programs 
used  and  of  critical  parts  of  the  experimental  statistical  analysis  and 
raw  data. 
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II  DATA  COLLECTION 


Previous  access-control  systems  using  hand  geometry  or  fingerprints 
have  used  devices  that  capture  data  from  physical  contact  with  or  proximity 
to  the  hand  surfaces  being  measured.  Such  contact  or  closeness  can  dis¬ 
tort  the  measurements ,  particularly  if  the  hand  can  move  or  exert  pressure. 

A  measurement  system  that  can  avoid  distorting  the  features  being  measured 
will  be  more  reliable,  which  in  turn  will  reduce  the  occurrence  of  incor¬ 
rect  authentication  decisions. 

Precise  orientation  of  the  hand  in  space  and  precise  control  of  finger 
and  thumb  spread  are  difficult  to  accomplish  without  some  mechanical  device, 
however.  To  compensate  for  the  variable  position  of  the  hand  from  trial  to 
trial  (because  there  was  no  device  to  hold  it  in  a  fixed  position)  it  was 
decided  to  concentrate  on  adjusting  the  data.  Precise  data  reduction  would 
allow  easier  and  more  natural  interaction  between  the  person  seeking  auth¬ 
entication  and  the  access-control  system. 

The  system  used  to  collect  data  is  schematically  illustrated  in  Fig¬ 
ure  1.  The  subject's  hand  is  positioned  near  the  light  diffuser;  the  open 
palm  faces  the  camera  and  the  fingers  are  comfortably  spread.  The  subject’s 
name  and  session  number  are  displayed  in  the  upper  quarter  of  the  imaged 
area  and  thus  become  an  integral  part  of  the  data  base. 

The  data  base  collected  contains  50  subjects,  25  males  and  25  females, 
ranging  in  age  from  the  early  20s  to  65,  in  height  from  approximately  5’ 
to  6' 3”,  and  in  weight  over  a  normally  observed  distribution.  The  total 
number  of  sessions  (five  trials  per  session)  was  1020,  approximately  evenly 
divided  between  males  and  females.  Of  these,  997  sessions  produced  4985 
high-quality  images  with  full  gray  scale.  In  the  remaining  sessions,  flash 
problems  of  failure  to  fire  or  synchronization  incurred;  the  resultant 
images  were  quite  useful  for  geometric  measurements  but  contain  no  gray¬ 
scale  features. 


FIVE  100W  SOFT-WHITE  BULBS 
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FIGURE  1  APPARATUS  USED  FOR  DATA  COLLECTION 
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Negative  images  have  been  cataloged  and  stored  in  three  volumes  that 
were  designed  and  manufactured  for  such  use.  Photo  development  followed 
standard  procedures  for  the  black-and-white  high-contrast  copy  film 
(Kodak  5069).  The  film  manufacturer  states  that  if  storage  conditions 
are  correct  (cool  and  dry)  the  archival  qualities  of  this  film  are  ex¬ 
cellent — it  will  retain  quality  for  more  than  30  years. 

Image  quality  is  also  excellent.  Correctly  focused  and  exposed  nega¬ 
tives  have  been  examined  with  a  400X  microscope,  and  minutiae  are  clearly 
visible  on  fingertips. 


Ill  DATA  PREPARATION 


A.  Data  Generation 

Two  methods  were  used  to  obtain  digital-image  data  to  be  used  for 
further  processing. 

•  An  SRI-built  Reticon  CCD  array  scanner  (512  elements  per  line) 
interfaced  to  an  LSI-11  was  used  early  in  the  project. 

•  A  commercially  available  Eikonix  digitizer  (1024-element  Reticon 
*  1500  scan  steps)  was  used  to  supply  final  evaluation  data. 

Data  from  the  first  system  were  generally  unsatisfactory  because  of 
a  combination  of  problems.  The  system  interface  included  a  hardware 
thresholder  for  the  analog  data  generated  by  the  Reticon  array.  The 
user  can  vary  the  threshold,  but  it  is  uniform  for  all  samples  coming 
from  the  CCD  array.  Such  an  arrangement  works  well  only  if  the  image  area 
is  uniformly  lighted,  the  individual  CCD  elements  are  reasonably  uniform  in 
sensitivity,  and  the  scanned  images  have  equivalent  niecom  densities. 
Unfortunately,  these  conditions  were  not  satisfied,  and  the  uniform  hard¬ 
ware  threshold  was  found  to  be  quite  unsatisfactory.  It  is  probably 
possible  to  build  a  system  that  utilizes  a  fixed-level  threshold  for  this 
type  of  application;  however,  such  a  device  with  a  sufficiently  uniform 

large  field  of  view  was  not  available,  and  the  uniform-threshold  technique 

* 

has  therefore  not  been  thoroughly  evaluate^ 

The  Eikonix  digitizer  is  an  excellent  instrument  for  purposes  of 
this  study.  The  digitized  images  generated  are  1024  *  1500-element  data 
each  with  8  bits  of  gray  scale  recorded  on  magnetic  tape.  The  real  reso¬ 
lution  of  the  data  when  the  film  negatives  are  sampled  with  this  spatial 
frequency  is  two  pixels  (picture  elements)  per  mm.  The  system  can 
automatically  compensate  for  lighting  nonuniformities  and  variable  CCD- 
element  sensitivities.  It  is  also  fairly  fast:  one  complete  cycle  of 
generating  and  recording  an  image  requires  about  one  minute.  All  of  the 
data  used  to  generate  the  results  appearing  in  this  report  were  collected 
through  the  Eikonix  system. 
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B.  Data  Reduction 


Image  data  on  magnetic  tapes  were  transported  to  a  PDP  11/40  computer 
system  in  SRI's  Bioengineering  Research  Center.  Each  image  was  reduced 
by  extracting  the  pertinent  512  *  512-element  field  and  by  thresholding 
to  a  binary  image.  Finally  the  binary  image  was  reduced  to  a  hand- 
perimeter  list  of  X-Y  coordinates  that  was  stored  on  a  DEC  RK  05  disk. 
Although  this  process  is  somewhat  involved,  it  does  provide  final  data 
that  are  quite  reliable  and  easy  to  work  with  in  subsequent  computations. 
(It  is  to  be  noted  that  the  final  data  set  generated  by  this  procedure 
does  not  include  any  Information  about  a  subject’s  hand  other  than  the 
perimeter  coordinates  in  a  512  x  512-element  field.) 

A  very  important  feature  of  the  data-reduction  process  is  the  method 
used  to  compute  the  perimeter  list  from  a  512  x  512-element  hand  image. 

The  method  involves  using  a  computer  program  (LINEBYLINE)  that  was  de¬ 
veloped  earlier  at  SRI  and  recently  rewritten  in  PASCAL  for  use  on  the 
PDP  11/40.  The  program  generates  the  perimeter  list  while  accessing  the 
image  data  one  line  at  a  time,  starting  at  one  edge  of  the  field.  This 
implies,  for  a  practical  system,  that  a  minimum  of  memory  and  data  buffer¬ 
ing  is  required;  consequently,  'there  is  a  significant  saving  of  money  and 
execution  time. 


IV  DATA  ANALYSIS 


The  basis  for  the  analysis  and  the  geometric  features  that  have 
been  extracted  is,  in  each  case,  a  list  of  the  X-Y  coordinates  of  the 
hand  perimeter.  Such  a  perimeter  is  described  with  a  real  resolution 
of  two  pixels  per  ram. 

All  measurements  are  based  on  landmarks  that  are  easily  identified 
on  the  perimeter;  the  tips  of  the  fingers  and  thumb  were  selected  for 
this  study.  These  landmarks  were  identified  by  calculating  a  coarse 
curvature  function  over  the  perimeter  and  finding  criterion  shifts  in 
the  angle  of  a  tangent  to  the  perimeter.  Tangents  were  estimated  by 
passing  an  18-point  window  over  the  perimeter  and  using  the  window  end 
points  to  define  a  vector  in  the  X-Y  image  plane.  The  window  was  moved 
in  nine-point  increments.  This  analysis  gives  a  rough  estimate  of  the 
tip  location. 

The  rough  estimate  is  refined  by  locating  four  additional  points 
on  the  finger,  two  found  25  ram  away  along  the  perimeter  on  each  side 
of  the  estimated  tip  and  two  more  found  35  mm  away  on  each  side.  The 
virtual  lines  connecting  each  pair  of  points  are  bisected,  and  the  inter¬ 
section  of  the  tip  and  the  line  defined  by  the  bisection  points  is  de¬ 
fined  to  be  the  actual  tip  location.  This  process  is  illustrated  in 
Figure  2.  After  a  precise  location  for  each  tip  is  established,  the 
following  features  can  be  extracted: 

•  Features  1  through  5  are  lengths  of  the  fingers  and  thumb, 
starting  with  the  little  finger.  Lengths  are  measured  from 
the  tip  to  the  depth  of  the  interdigital  space  as  illustrated 
in  Figure  3. 

•  Features  6  through  10  are  widths  of  the  fingers  and  thumb. 

A  width  is  the  distance  between  two  points  on  the  perimeter 
of  a  finger;  the  points  are  45  mm  down  the  perimeter  on  each 
side  of  the  tips  of  the  little  finger  and  the  thumb,  and  60  mm 
on  each  side  of  the  other  t  jertips.  Width  measurement 
points  are  illustrated  in  Figure  4.  These  features  are 
ordered  as  were  the  first  five. 


10mm 


FIGURE  2  FINGERTIP  LOCATION  MEASUREMENT 


FIGURE  3  MEASUREMENT  POINTS  FOR  LENGTHS 
OF  FINGERS  ANO  THUMB 
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FIGURE  4  FINGER  AND  THUMB  WIDTH  MEASUREMENT 
POINTS 


•  Features  11  through  15  are  ratios  of  the  measurements  of  finger 
lengths  to  widths. 

•  Feature  16  is  the  hand  width  measured  between  two  points  on  the 
palm.  One  point  is  9  cm  on  the  perimeter  from  the  tip  of  the 
little  finger  toward  the  wrist  and  the  other  is  11  cm  on  the 
perimeter  from  the  tip  of  the  index  finger  toward  the  thumb. 
These  points  are  labeled  2  and  11,  respectively,  in  Figure  4. 

•  Features  17,  18,  and  19  are  the  hand  perimeter,  hand  area,  and 
the  ratio  of  the  squared-perimeter  length  to  the  area.  The 
perimeter  and  area  are  measured  from  point  1  to  point  14  in 
Figure  4.  Point  1  is  13  cm  on  the  perimeter  from  the  tip  of 
the  little  finger,  and  point  14  is  10  cm  on  the  perimeter  from 
the  tip  of  the  thumb. 

•  Ten  features,  20  through  29,  are  paired  values  corresponding 
to  the  curvature  of  each  fingertip  and  the  distance  over  which 
the  tip  shape  can  be  represented  by  a  circular  arc  with  less 
than  a  criterion  value  of  rms  (root  mean  square)  error.  The 
process  of  extracting  a  typical  feature  pair  is  illustrated  in 
Figure  5.  A  tangent  to  the  perimeter,  such  as  T,  has  an  angle, 
a,  as  shown.  If  a  is  plotted  as  a  function  of  distance  along 
the  perimeter  defining  a  finger,  then  a  plot  such  as  the  one 
shown  in  Figure  5  results.  The  slope  of  the  plot  between  P-^ 
and  P2  is  related  to  the  curvature  of  the  tip.  The  slope  is 
recorded  as  the  curvature.  A  correlation  calculation  over  the 
region  containing  this  characteristic  region  yields  a  high  value 
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PERIMETER 


FIGURE  5  CURVATURE  AND  DISTANCE  MEASURES 


until  the  boundary  of  the  region  extends  into  the  flat  segments 
on  either  side  of  the  sloping  section.  Pj  and  P2  are  the  boundary 
points  when  the  correlation  coefficient  (using  a  straight-line 
fit  to  the  sloping  region)  falls  to  0.98.  The  distance  between 
P-^  and  P2  is  recorded  as  the  perimeter  length  over  which  the 
fingertip  can  be  approximated  well  with  a  circular  arc. 

In  summary,  hand  geometry  is  expressed  as  a  set  of  29  features  that 
are  measured  by  locating  precisely  five  landmarks  on  the  perimeter  of  the 
hand — the  tips  of  the  fingers  and  the  thumb.  The  feature  set  contains 
16  measures  of  lengths  and  widths  and  ratios  of  lengths  and  widths,  and 
three  measures  of  area  or  area  and  perimeter  length.  Ten  measures  of 
shape  (curvature)  are  included. 


V  CLASSIFICATION  RESULTS 


Thirty  subjects  (15  male,  15  female)  were  selected  randomly  from 
the  subject  population  and  278  images  taken  from  the  selected  subjects 
were  digitized.  The  data  were  converted  to  perimeter  lists  and  sets  of 
29  features  were  extracted.  The  data  were  submitted  to  a  standard  sta¬ 
tistical  analysis  program*  to  determine  the  power  of  the  extracted 
features  in  discriminating  among  the  30  subjects.  Linear  discriminant 
analysis  was  performed  and  all  278  feature  sets  were  correctly  classi¬ 
fied  for  all  subjects.  The  detailed  results  of  this  analysis,  as  well 
as  the  feature  data  set,  are  submitted  as  a  separate  package.  Annotation 
for  this  package  is  included  as  Appendix  C. 

Of  critical  importance  are  indicators  of  value  and  reliability  of 
the  various  features  selected.  One  such  indicator  is  the  univariate 
F-ratio,  a  number  that  expresses  the  ratio  of  be tween- groups  variation 
of  a  given  measure  and  the  within-groups  variation  of  that  same  measure. 
Large  values  of  F  indicate  that  a  measure  will  be  relatively  useful  in 
separating  groups  (or  samples  from  different  subjects)  and  small  values 
predict  that  a  measure  will  do  poorly.  Figure  6  summarizes  the  F-ratios 
of  the  selected  29  features  over  the  30-subject,  278-sample  data  base. 

All  of  the  features  appear  to  be  useful  and  reliable,  especially 
the  more  global  measures  of  hand  width,  perimeter,  and  area.  Some  mea¬ 
sures,  notably  the  curvature  and  distance  measures  associated  with  the 
thumb  are  marginal  as  presently  calculated.  They  could  be  made  more 
useful  by  improving  the  algorithm  that  deals  with  the  curvature  of  the 
perimeter  in  the  vicinity  of  the  thumb  and  by  locating  the  thumb  tip 
more  precisely. 
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.0709 

113.0063 

FEATURU 

.3051 

19.9798 

FEATUR16 

.0282 

295.1276 

FEATUR17 

.0185 

953.6191 

FEATURIB 

.0152 

559.9787 

FEATURU 

.0950 

181.5809 

FEATUR20 

.9010 

12.7763 

FEATUR21 

.1229 

61.0282 

FEATUR22 

.1865 

37.3090 

FEATUR23 

.2269 

29.2203 

FEATUR24 

.1555 

96.9599 

FEATUR25 

.1325 

55.9899 

FEATUR26 

.9777 

9.3999 

FEATUR27 

.2098 

33.1988 

FEATUR2B 

.  7087 

3.5192 

FEATUR29 

•  5388 

7.3202 

FIGURE  6  ANALYSIS  OF  HAND  GEOMETRY  FEATURES 


In  general,  it  appears  that  the  measures  used  are  being  generated 
very  reliably,  which  results  in  good  (high)  F-ratios.  One  difficulty 
in  previous  systems  that  used  hand  geometry  for  authentication  was  un¬ 
reliable  measures  of  certain  features,  with  correspondingly  small 
F-ratios  and  poor  discriminating  power. 

Figure  7  gives  more  information  about  the  quality  of  the  feature 
set.  In  this  table,  the  within-groups  correlation  matrix  is  listed. 

The  correlation  matrix  is  used  to  express  the  degree  to  which  pairs  of 
features  are  correlated  or  dependent.  Circled  are  the  two  largest  mag¬ 
nitudes  of  correlation  values,  corresponding  to  the  pairs 

•  Second  finger  length  and  the  ratio  of  its  length  and  width 

•  Thumb  width  and  the  ratio  of  its  length  and  width. 

It  is  to  be  expected  that  ratios  of  features  and  the  features  themselves 
will  tend  to  have  relatively  large  correlation  magnitudes,  and  the  data 
in  Figure  7  substantiate  this  expectation.  However,  the  SPSS  results 
show  that  all  features  contribute  positively  in  the  discriminant  analysis. 

In  summary,  we  have  been  able  to  improve  on  previous  attempts  to  use 
features  of  the  hand  for  access  control  In  the  areas  of  accuracy  and  user 
comfort.  This  has  primarily  been  the  result  of  removing  the  requirements 
on  both  the  subject  and  the  mechanical  device  that  produced  either  a 
distortion  of  the  data  or  discomfort  for  the  individual.  Instead,  we 
have  required  the  software  to  do  all  the  registration  and  measurements. 

The  result  is  a  fast  and  reliable  method  that  demands  only  a  minimum  of 
effort  from  the  individual  requesting  access. 

If  this  technique  were  to  be  packaged  into  a  commercial  product, 
we  estimate  that  the  throughput  rate  could  be  less  than  two  seconds. 

The  rate-limiting  step  in  the  process  is  data  acquisition  from  the 
scanner  through  an  interface  and  into  a  computer.  The  throughput  rate 
of  a  practical  interface  could  be  in  excess  of  275K  bytes  per  second. 
Because  each  image  contains  256K  bytes,  the  transfer  would  take  about 
one  second.  This  data  rate  is  well  within  the  state  of  the  art  of 
image  production  using  a  scanning  linear-array  system.  Because  the  data 
can  be  analyzed  a  line  at  a  time,  and  therefore  analyzed  as  it  is 
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FIGURE  7  ANALYSIS  OF  HAND  GEOMETRY  FEATURES  USING  WITHIN-GROUPS  CORRELATION  MATRIX 
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produced,  we  do  not  believe  that  the  data  processing  would  be  limiting. 
However,  if  line  processing  takes  longer  than  scanning,  multiple  proces¬ 
sors  could  be  used  with  parallel  multiple  lines  of  processing  working 
concurrently  to  process  more  than  one  line  at  a  time.  The  need  for 
multiple  processing  has  not  been  assessed,  since  all  work  has  been  per¬ 
formed  on  a  multi-user  system  with  relatively  slow  development- type  pro¬ 
grams.  However,  because  our  software  could  be  supported  by  multiple 
processors  working  in  parallel,  process  time  should  not  be  rate-limiting. 

Since  the  potential  of  this  technique  has  been  shown  to  be  good, 
the  next  task  is  to  refine  the  technique*  Refinements  can  be  made  in 
all  areas,  including  fine-tuning  feature  extraction,  adding  new  features 
and  removing  the  less  significant  ones,  and  developing  methods  for  making 
authentication  decisions  based  on  these  data.  In  addition,  studies 
could  be  made  of  methods  for  optimizing  the  software  for  minimum  through¬ 
put  time . 
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VI  COMPUTER  PROGRAM  SUMMARIES 


The  task  of  image  analysis  has  been  broken  down  into  four  discrete 

jobs: 

•  Determine  the  correct  gray  level  for  image  thresholding  using 
the  FORTRAN  program  HIST. 

•  Threshold  the  image  to  produce  a  binary  image  using  the  FORTRAN 
program  THRESH. 

•  Extract  a  list  of  hand-perimeter  points  from  the  binary  image 
using  the  PASCAL  program  LINEBYLINE. 

•  Extract  the  feature  values  from  the  perimeter  list  using  the 
FORTRAN  program  SCAN. 

Each  of  these  jobs  is  handled  by  a  single  independent  program.  A  brief 
description  of  each  program  follows,  and  listings  of  each  program  appear 
in  Appendix  A.  The  program  LINEBYLINE  is  quite  complex,  however,  and 
Appendix  B  is  included  to  describe  its  operation  in  detail. 

A.  HIST 

This  FORTRAN  program  computes  a  histogram  of  the  frequency  of  occur¬ 
rence  of  gray  levels  in  an  8-bit  gray-scale  image.  The  program  is  capable 
of  handling  a  512  x  512  pixel  image,  which  is  the  size  used  to  represent 
the  hand  data.  The  image  format  is  two  pixels  per  16-bit  word.  The  record 
size  on  a  DEC  RK  05  disk  is  256  16-bit  words.  The  output  file  lists  the 
frequency  count  of  each  gray  level  (0-255).  A  plot  of  the  output  can  also 
be  produced  on  a  TEKTRONIX  4000-series  terminal,  using  the  PLOT  10 
(TEKTRONIX)  software  interface  installed  in  the  Bioengineering  Research 
Center's  PDP  11/40  computer.  The  output  of  this  program  clearly  'identi¬ 
fies  the  sets  of  gray  levels  that  represent  either  the  background  or 
the  hand;  this  information  is  used  to  pick  a  discrete  gray  level  between 
the  two  sets  of  levels  that  will  best  separate  them.  This  level  is  used 
by  THRESH  to  produce  a  binary  image. 
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B.  THRESH 


This  FORTRAN  program  thresholds  an  8-bit  gray-level  image  to  produce 
a  binary  image.  Thresholding  is  performed  on  a  pixel-by-pixel  basis. 

The  input  file  contains  two  8-bit  pixels  per  word  with  a  total  image  size 
of  262,144  pixels.  The  output  file  contains  16  one-bit  pixels  per  word. 
Both  input  and  output  files  have  256  words  per  record.  The  discrete 
threshold  value  is  entered  at  run  time. 

C.  LINEBYLINE 

This  PASCAL  program  produces  a  list  of  perimeter  points  that  bound 
"blobs"  of  color  opposite  to  that  of  the  background  in  a  binary  image. 
Eight-point  connectivity  is  checked  for  each  point  in  the  image.  Using 
the  large  area  as  a  criterion,  the  hand  blob  is  selected  from  others, 
and  only  the  hand  perimeter  is  saved.  This  list  of  hand-perimeter  points 
is  output  as  a  file  of  ordered  points  running  around  the  perimeter  of  the 
hand.  Only  one  pass  is  made  through  the  image,  and  only  one  line  is 
accessed  at  a  time.  The  program  therefore  requires  only  slightly  more 
dynamic  memory  than  is  required  to  store  perimeter  points.  A  more  exten¬ 
sive  description  of  this  program  is  included  as  Appendix  B. 

D.  SCAN 

This  FORTRAN  program  processes  the  list  of  perimeter  points  produced 
by  LINEBYLINE  and  outputs  numeric  values  that  describe  geometric  features 
of  the  hand,  as  discussed  earlier  in  this  report.  Included  below  are 
comments  intended  to  help  the  reader  interpret  the  program  listing  in 
Appendix  A. 

•  The  preferential  landmarks  on  a  perimeter  are  the  fingertips. 
These  are  located  by  the  angle  algorithm  on  page  A-20, 
starting  with  the  comment  lines.  Landmarks  are  found  by 
moving  a  window  of  size  IDIS  (18  points)  along  the  perimeter 
and  looking  for  an  angle  criterion  shift  of  90°  in  the  tangent 
to  the  perimeter. 

•  Starting  at  line  174  on  page  A-5,  calculations  of  finger  length 
are  made.  These  require  the  precise  tip  locations  calculated 
earlier  on  page  A-4,  starting  with  the  comment  lines  in  the 
middle  of  the  page. 


•  Finger  widths  are  calculated  with  the  code  on  page  A-21. 

•  Ratios  of  finger  lengths  to  widths  are  calculated  on  page  A-22. 

•  Curvature  measures  are  calculated  with  the  code  in  the  middle 
of  page  A-22,  using  subroutines  TANAG  and  LINFIT  that  process 
angle-function  data. 

•  Hand  width,  perimeter,  and  hand  area  are  calculated  at  the  bottom 
of  page  A-22,  using  the  functions  IDELTS  and  DSTNCE  and  the  sub¬ 
routine  ARPRIM. 


21 


Appendix  A 

COMPUTER  PROGRAM  LISTINGS 


nnnn  oooo o  ooooooooc*  v*  oo  oo  09  o  00  ooonoonooonn 


PROGRAM  HIST 


WITTEN  BY  GREG  MYERS 

THIS  ROUTINE  COMPUTES  A  HISTOGRAM  OP  AN  8-BIT  GREY-LEVEL  IMAGE.  IT 
HANDLES  IMAGES  THAT  ARE  256,  512,  OR  1624  POINTS/LINE.  THE  IMAGE  FORMAT  IS 
2  PIXELS  PER  WORD,  AND  THE  RECORDS l ZE  ON  DISK  IS  256  16-BIT  WORDS. 

THE  OUTPUT  FILE  LISTS  THE  FREQUENCY  COUNT  OF  EACH  GREY  LEVEL  (0-2SS) . 

NOTE:  FOR  256  X  256  FILES,  THERE  IS  ONLY  ONE  PIXEL  PER  WORD! 


LOGICAL*  l  NULL ,  NAM  IN  (30),  NAMOUTO0) ,  A<5I2>,  ITITLE<72> 
l NT EGER* 4  IH<2S6) ,  IHAKEA 

REAL  G(2S6),  H(2S6),  XAXISU8),  YAX1S12)  ,  TITLE! 18) 
INTEGER  STARTL,  STARTP ,  ENDL,  ENDP 
EQUIVALENCE  < IT ITLE.TITLE) 

DATA  XAXISl  1  >/ ’GREY '/,XAXIS(2)/ '  LEVV.XAXISOi/'EL  V, 

I  YAXISI  1)/‘C0UNV,YAXIS(2)/%T  '/ 


LUNIN  -  3 
Luiioirr  •  4 
NULL  •  0 

DO  I  I  *  4,18 
TITLE! I)  -  4H 
I  XAXIS(I)  •  4H 

*— GET  DATA  FROM  TERMINAL 


WRITE  (5  10) 

10  FORMAT < •  FOR  INPUT  FILE:  TYPE  0LINES,  # POINTS/LINE,  AND  FILENAME' > 
READ  (5,15)  NL.  NP,  NCHIN ,  (NAM1N< I ) , I- I .NCHIN) 
is  format<215.q,30Ai> 

NAMlN(NCHIN^l)  •  NULL 

OPEN(UNIT-LUN IN .NAME-NAN IN, TYPE* ’OLD*  .ACCESS* ’DIRECT’  .READONLY, 

1  RECORDS  IZE*  128,  ERR -901) 


WRITE  (5,20) 

20  FORMAT < ’  FOR  HISTOGRAM  AREA:  TYPE  STARTING  LINE,  STARTING  POINT,  I 
1LINES,  #POINTS/LlNEM 

READ  (5.25)  STARTL,  STARTP,  ML.  MP 
25  FORMAT (415) 

ENDL  •  STARTL  ♦ML  *  I 
ENDP  •  STARTP  ♦  MP  -  l 


50 


55 

40 

45 


1 

2 


I 


WRITE  (S  30) 

FORMAT  < '  *TYPE  *’P*  IF  YOU  WANT  A  HISTOGRAM  OUTPUT  FILE  ONLY,'/ 
'  "2"  IF  YOU  WANT  A  PLOT  ONLY,*/ 

»’3"  IF  YOU  WANT  BOTH’) 

READ  (5,35)  ISVP 
FORMAT  (ID 

IF  <ISWP  .EO.  I  OR.  ISWP  .EQ.  3)  WRITE  (5,40) 

FORMAT <  *  TYPE  HISTOGRAM  OUTPUT  FILE’) 

IF  (ISWP  .EQ.  1  OR.  ISVP  .EQ.  3) 

READ  (5,45)  NCXIOUT, <NAMOUT( I ) , I-l .NCHOUT) 

FORMAT  (Q,30AU 


NAMOUT (NCHOUT- 1)  •  NULL 

IF  (ISWP  .EQ.  2  .OR.  ISWP  .EQ.  3)  VRITE<5,56) 

8  FORMAT ( ’  TYPE  0  OF  GREY  LEVELS  DESIRED  IN  THE  PLOT  (13)') 

IF  (ISWP  .EQ.  2  .OR.  ISWP  .EQ.  3)  READ(5,55)  NGREYL 

5  FORMAT (13) 


—  INITIALIZE  HISTOGRAM  ARRAY 


106 


TIN  -  SECNDS(6.) 
DO  106  I  -  1,256 

ih< i )  •  a 


IF  (NP  .0) 

IF  (NP  .05 
IF  (NP  .0) 
GO  TO  962 


1624)  GO  TO  260 
512)  00  TO  360 

256)  00  TO  466 


— FOR  1624  POINTS/LINE  (2  RECORDS  PER  LINE) 


3M  DO  210  JRBC  -  <SMHLM>«2M1<RNBlrl>«2»laft2 

D  IF  (NOOURGC.20)  .BO.  1)  WRITE  (5,212)  JRBC 

D  212  FORMAT*  *  JRBC  •  9 ,14) 

IF  (STARTF  .GT.  512)  00  TO  220 

READ  <UmnrjREC,ERR-903>  (A< I) f  1*1 .512) 

K1  -  STARTF 

12  •  MIN0I ENDP, 5)2) 

CALL  HISTL(A,Kt ,K2, !H) 

220  CONTINUE 

IF  (ENDP  .LE.  512)  00  TO  210 

READ  ‘LUNIN 1 JRBC*  1  ,  ERR  *904)  (Ad) , !•! ,S|2> 

K 1  -  MAX0(STARTP,S13>  -  512 

K2  -  ENDP  -  512 

CALL  HISTL(A,KI ,K2, IH) 

210  CONTINUE 

oo  ro  see 
c 
c 

C - FOR  512  POINTS/LINE  (ONE  LINE  FER  RECORD) 

C 

300  DO  310  JRBC  •  STARTL, ENDL 

READ  (LUNIN* JREC.ERR-903)  <A( I ) . I ■ 1 ,512) 

CALL  HISTLU, STARTF, ENDP.  IH> 

310  CONTINUE 

00  TO  500 
C 

C  _ 

C - FOR  256  POINTS/LINE  (2  LINES/RBOORD) 

C 

400  DO  410  JREC  -  STARTL , FNDL 

READ  (LUNIN* JREC.ERR-903)  <A( 1 ) , !• 1 ,512) 

C— COMPRESS  DATA  IN  ARRAY  A  FROM  DEANZA  FORMAT 
DO  411  I  -  1,256 
411  A<l)  -  A*2*|- 1 > 

CALL  HINTLiA.STARTP.ENDP, IH> 

410  CONTINUE 
C 
C 
C 
C 

C - PRODUCE  OUTPUT  FILE  IN  A  FORMAT  ACCEPTABLE  TO  THE  PROGRAM  PLOT  DAT  A 

C 

500  CONTINUE 

TINE  •  SBCH0S<T1H> 

WRITE  (S,50|)  TIME 

501  FORMAT <  ’  TIME  *  \H0.3,*  SECONDS1) 

CLOSE  (UNIT-LUNIN) 

DO  565  1  •  1,256 
Gd)  -  1-1 
H<i>  •  m< i) 

505  CONTINUE 
C 

IF  ( ISRP  .NE.  1  .AND.  ISVF  .NE.  3)  00  TO  600  _ 

OPEN  <  UN  IT-LITNOUT ,  NAME*  MAHOUT ,  TYPE-  *  NEW  \ACCESS-  ’SEQUENT  I AL  * , 

|  FORM* 1  FORMATTED* , ERR -905. CARR  1  ACBOONTROL-xIsT •  > 

C 

WRITE  (4,510)  (NAM IN ( I > . I- 1 ,27> .  STARTL, STARTF .ML, NP 
516  FORMAT  (7X*  HISTOGRAM  OF  \27Al,4fs> 

N256  »  256 

WRITE  (4,520)  N2S6,  <0( !> ,H( I) • 1*1,256) 

520  FORMAT ( I4,/,2(E15.8,2X) > 

I H AREA  -  0 
DO  530  I  *  1,256 
530  IHAKEA  *  IHAREA  ♦  111(1) 

WRITE  (4,540)  IHAREA 
540  FORMAT!/.*  HISTOGRAM  AREA  -  *,17) 

CLOSE  (UN1T»LUN0UT) 

C 

C 

C - PLOT  THE  DATA  ON  THE  TEKTRONIX 

C 

600  CONTINUE 

IF  (ISWF  LT.  2)  GO  TO  909 

ENCODE  (67,510,  ITITLE,  ERR-9)  (NAMINd ) ,  !•  1 ,27)  , STARTL, STARTF, ML, I 
CALL  PLOT  <G,H,NGREYL,XAX  IS,  YAXIS,  TITLE) 


WRITE  (5,911)  (NAMINd).!- 1,1 NCHIN) 

FORMAT ( *  HIST:  ERROR  OPENING  INPUT  FILE  *,30AD 
00  TO  999 
WRITE  (5,912) 

FORMAT ( *  HIST  ERROR:  #  POINTS/LINE  256,  512,  OR  1024  ONLY t ’ ) 
CLOSE  (UNIT-LUNIN) 

00  TO  999 
WRITE  (5,913) 

FORMAT ( ’  HIST  ERROR:  ERROR  903  READING  INPUT  FILE*) 

CLOSE  (UNIT-LUNIN) 

00  TO  999 


oonooooooooo  oooonoooooononoooo 


964 

914 


90S 

9tS 

9 

916 


C 

C 


106 

C 

C 

999 

C 

C 


WRITE  (5,914) 

FORMAT ( *  HIST  ERROR:  ERROR  964  READING  INPUT  PILE') 
CLOSE  (UNIT-LUNIN) 

00  TO  999 

WRITE  (5,9151  <NAMOVT<I).I-l,NCIIOtnr) 

FORMAT ( ’  HIST  :  ERROR  CREATING  OUTPUT  PILE  \30A1> 
GO  TO  999 

WRITE  (5.916)  (ITITLE(l), 1-1,72) 

FORMAT ( 72A 1 ) 

PAUSE  916 
00  TO  999 

END 

SUBROUTINE  HISTL<A,KI ,K2, !fl> 

UDGICAL*1  AC512) 

INTEGER*4  1H<256> 


DO  100  I  •  K1.K2 
INT  •  Ad) 

IF  (INT  , LX.  0)  INT  •  INT  ♦  256 
1NDX  -  INT  ♦  l 

IF  (IH(INDX)  .EQ.  32767)  00  TO  106 

IH(INDX)  •  IH(INDX)  *  1 

CONTINUE 


RETURN 

END 


PROGRAM  THRESH 


WRITTEN  BY  ORB!  MYERS 


THIS  PROGRAM  THRESHOLDS  AN  8-BIT  GREY-LEVEL  IMAGE  AND  PRODUCES  A  BINARY 
<1  B|T  PER  PIXRL>  OUTPUT  IMAGE.  THE  INPUT  FILE  CONTAINS  l  8-BIT  PIXEL 
PER  WORD  IF  THE'  IMAGE  IS  IN  THE  256  X  256  DE  ANZA  FORMAT.  OTHERVIZE, 
THE  IMAGE  CONTAINS  2  PIXELS  PER  WORD.  THE  OUTPUT  FILE  CONTAINS 
16  1-BIT  PIXELS  PER  WORD.  BOTH  THE  INPUT  AND  OUTPUT  FILE?  HAVE 
256  WORDS  PER  RECORD. 


MODIFIED  9/24/86  TO  ACCEPT  COMMAND  STRING  INPUT 
COMMAND  SYNTAX  IS: 

I OUTF ! L ) ■ I NF 1 L/LN : * : #/TH : # 

ROB  HAMMON 

LOGICAL* 1  NULL,LFLG,NANFIL<34) , A(8192) ,B(512) 
INTEGERS  I  BYTE 
LOGICAL* I  BYTE 

INTEGER  1FLG,NL,NP ,  1THR.EXSTAT 
INTEGER  SVTAB(9),THVAL(S),LNVAL<5) 

REAL  PROMPT 

DATA  PROMPT/’ THBV,  IEXT/3REXT/.EXSTAT/1/ 

DATA  NL/-1/.NP/-1/.ITHR/-1/ 

EQUIVALENCE  (BYTE, 1 BYTE) 

DATA  LUNIN/3/,  LUNOUT/4/,  NULL/0/ 


— GET  DATA  FROM  TERMINAL 


16 

16 


WRITE  (5,16) 

FORMAT ( '  FOR  INPUT  FILE:  TYPE  ALINES,  • POINTS/LINE,  AND  FILENAME’) 
READ  (5.15)  NL,  HP,  NCH1N,  (NAH1N( 1) , I- 1 ,NCHIN) 

FORMAT (ilS,Q,36AI) 

NANINOfCHIN* 1)  •  NULL 


Iaitialilo  GETCMD  and  twitoh  datcriptar  tttblu 

CALL  ASSIGN  <6,’TI:\3>  fAttiga  LUN  for  GETCMD 

CALL  INICND  (6, ! EXT, PROMPT)  !S«t  LUN  and  fila  ait, prompt 

CALL  CS1SW  ISVTAB( I > , ’LN* ,2, , LNVAL)  fDaolara  AM  twi loh 

CALL  CS1SW  (SVTAB(S) , ’TH* , I , ,THVAL)  fDaolara  /TH  twitoh 

CALL  CS1SV  (LNVAL(1),’D’,NL>  fDaolara  Mrtt  *al  for  /LN:i 

CALL  CSISV  (LNVAL(3),’D’,NP)  tSoeoad  volte 

CALL  CSISV  (THVAUI),’D\ITHR>  fDaolara  ral  for  /TH:# 


C  0#t  >0— m4  string  pr«o«ia  it 


16 


♦Error 

tStl  sicotnsf  n 1 

♦Coaprts*  string 
♦Error 

♦Eqtnl  sign  it  string 
♦No  tqoal  sign 
*S#t  deftnlt  tit 
♦  St  t  for  ohscli 


♦Pnrst  string 


♦Error 


CALL  GETOtD  (B, , , .  ,LFLG> 

IF  <LFLO  .NE.  6)  00  TO  1906 
EXSTAT  •  1 
CALL  CSI 1 (B,  t , 1FL0) 

IFIIFLG  .GT.  0)  00  TO  1926 
IFUFLG  .EX).  0>  LFLO  •  M* 

IFUFLG  .LT.  0>  LFLO  •  *0* 

CALL  CS1DEF(LFLG, , , , *B8  ’) 

NL  •  -1 
NP  •  - 1 

inm  •  -1 

CALL  GSI2(LFLG,NANPIL«iSVTAB) 

IF(NAW7tL<U  BO.  O)  GOTO  2949  fError 

IFUNL  .EO.  -I)  OH,  (NP  .BQ.  -U>  GOTO  1930 
IFtITHR  .BQ.  -1)  GOTO  1030  fError 

OPEN  <  UN  I T-LUN  I N ,  NAME»NAMF  1 L ,  TYPE*  *  OLD  *  .ACCESS*  'DIRECT' ,  READONLY , 
RECORDS I ZE* 120,ERI*90l ) 


TOTPIX  ■  FLOAT (NL)  •  HP 
NRECIN  -  TOTPIX  /  SI2 
NBLOUT  •  TOTPIX  /  4996 
L1NC  •  1 

rror  •  9 

IF  (NL  .HE.  256  .OR.  NP 
NRECIN  •  NRECIN  •  2 
LI NC  *  2 
KTOT  •  16 
CONTINUE 


HE.  256)  00  TO  10 


19 

C 

C  Gst  inforastion  tnd  opts  otiptt  ft  It 
IF  (LFLG  .BQ.  *1*1  00  TO  59 
NANFIL(2S>  •  *B* 

NABF1L(26>  -  M* 

NANFIL(27)  •  ’N* 

GO  TO  69 

56  CALL  CSIDEF( *0’ , , ,NANF1L( IS) , 'BIN* ) 

CALL  CSI2( *0’ .NANFIL) 

IP(NANFILU)  .EO.  0)  00  TO  1049 
WRITE  (5  20) 

FORMAT ( '  TYPE  THRESHOLD,  AND  NARE  FOR  BINARY  OUTPUT  FILE*) 
READ  (5,25)  ITHR,NCHOUT,(NANOUT(|) , Nl.NCHOUT) 

FORNAT ( 14 ,Q, 39A I ) 

N ABOUT (NCHOUT* I )  -  NULL 


♦Eqnel  sign  is  string 
♦Stt  dtftolt  tttaasion 


♦Set  dtftolt  sit 
♦Ptrst  string 


C 

C  29 
C 

C  25 
C 

66 

C 

C 

c 

c 


l 


OPEN  ( UN  I T *LUNOUT ,  NANE-NAMF IL ,  TYPE-  *  NEW  * ,  ACCESS*  *  D I RBCT  * , 
RECORDS IZE- 128, ERR *902, 1NIT ULS IZE* NBLOUT) 


TIB  -  SBCNDSC9.) 

DO  109  J  *  1 .NRECIN, KTOT 


C - READ  IN  4096  0-BIT  PIXELS 

C 


299 


C 

C 


C 

C- 

c 


DO  209  K  •  1.ICT0T 

KOFF  •  (K-l)*512 
READ  (LUNIN’ J»K-l) 

CONTINUE 


LOFF  ■  I 

DO  390  L  -  1.512 
I BYTE  -  0 


(A(KOFFtl) ,1*1,512) 


-PERFORM  THE  THRESHOLD  IN  GROUPS  OF  0  PIXELS 


DO  400  I  «  1,0 

l NT  *  A(LOFF) 

LOFF  *  LOFF  ♦  LINO 


♦  256 


399 


eon  tints 


IF  (INT  .LT.  9)  INT  -  INT 
ibytt  •  «shf t< ib|tn,-l) 
if  <int  .gt.  ithr)  ibyta  ■  tor ( Ibytt ,  120) 

con t its# 

b ( 1 )  -  by tn 


t - wrl tn  ont  ent  rtoord  of  dots  (4996  btnory  pistls) 

writs  Unto*  t  *  j/X  tott  1)  (b  ( I ) ,  i«  1 ,512) 

190  ootlisto 


CLOSE  (UNIT  -  LUNIN) 
CLOSE  (UNIT  -  LUNOUT) 
TINE  ■  SECNDS(TIN) 

WRITE  (5,991)  TINE 

forratc  tire  •  \Fie.3> 


A-5 


999  00  TO  10 


tGnt  ■•it  eoaaiid 


901  writ*  (5,911)  (NANFILI i) , i»l .34) 

911  fopiitC  •rror  op«aii|  lipat  ril#  ’,30al> 
to  to  999 

902  writ*  (5,912)  (NAMFILU >  .  i- 1 ,34) 

9)2  forantC  #rror  ornatiat  oitpit  filn  ’,38*1) 
to  to  999 

C 

1000  IF(LFLO  .EQ.  366)  CALL  EXST(EXSTAT)  fExit  with  ttatm 

1020  TYPE  1021 

1021  FORMAT  <’  ###  eouud  ayatax  orror  •••*> 

00  TO  9000 

1030  TYPE  1031 

1031  FORMAT  (’  •••  orror  in  inpit  fitospno  •••*) 

GO  TO  9600 

1040  TYPE  1041 

1041  FORMAT  <’  •••  Error  in  oitpit  filnapne  ••••) 

00  TO  9000 

C 

9000  EXSTAT  ■  4  tS*t  for  aa?#r  nrror 

—  CALL  RESCMD  fR«s«t  cad  input  to  top  lovol 

00  TO  10  10*1  anxt  ooanind 


nod 


PROGRAM  LINEBYLINE(TTY);  (•NSP«> 


(•  VERSION  46C  10/16/80 

THIS  VERSION  PRODUCES  A  PERIMETER  LIST  OF  ONE  LARGE  BLOB,  WHICH  IS  SUPPOSED 
TO  BE  A  HAND.  FOR  ROB  HANNON  AND  THE  PALNPRIMT  STUDY. 

WRITTEN  BY  GREG  MYERS 

SRI  INTERNATIONAL,  J-3006 
333  RAVENSVOOD  AVENUE 
MENLO  PARK,  CA.  94025 


THIS  ALGORITHM  IS  TAKEN  FROM  A  REPORT  BY  GERRY  AGIN  AT  SRI  ENTITLED 
IMAGE  PROCESSING  ALGORITHMS  FOR  INDUSTRIAL  VISION.  THIS  ALGORITHM 
IS  REFERRED  TO  AS  CONNECTIVITY  ANALYSIS" .  IT  SEGMENTS  A  BINARY 
IMAGE  INTO  "BLOBS"  (CONNECTED  AREAS)  OF  THE  SAME  COLOR  (GREY  LEVEL). 

ONLY  ONE  PASS  IS  MADE  THROUGH  THE  IMAGE,  AND  ONLY  ONE  LINE  IS  ACCESSED 
AT  A  TIME  (HENCE,  THE  PROGRAM  NAME  ’LINE-BY-LINE’ > .  FEATURES  OF  EACH 
BLOB  ARE  COMPUTED,  SUCH  AS  ITS  AREA,  CENTER  OF  GRAVITY,  BOUNDING  RECTANGLE, 
AND  A  LIST  OF  ITS  PERIMETER  POINTS. 


THE  LETTERS  NSP*  WITHIN  A  COMMENT  MEAN  ’NON-STANDARD  PASCAL*. 

THE  LINES  NARKED  ’NSP’  NAY  REQUIRE  MODIFICATION. 

THE  LETTERS  ’NO'  WITHIN  A  COMMENT  MEAN  ’MACHINE-DEPENDENT'. 

THESE  LINES  NAY  ALSO  REQUIRE  MODIFICATION. 

BECAUSE  THE  COMMAND  ’DISPOSE’  IS  NOT  IN  THE  SWEDISH"  VERSION  OP  PASCAL 
AT  SRI,  DISPOSE’  STATEMENTS  ARE  ENCLOSED  IN  COMMENT  STATEMENTS  AND 
ARE  IGNORED. 

READ  AND  WRITE  STATEMENTS  FROM  THE  USER’S  TERMINAL  DO  NOT  REQUIRE 
FILE  SPECIFICATION  IN  THIS  VERSION  OF  PASCAL.  THE  FILE  ’TTY*  IS 
ASSUMED. 

EDLN  IS  AT  THE  BEGINNING  OF  AN  INPUT  LINE  WHEN  READING  FROM  THE  TTY. 

*) 


CONST  NPIXELSPERLINE  •  512; 
NUNES  •  512; 

NB1TSPERC1IAR  *  8: 
NCIlARSPFRL I NE  -  64:  (• 

NCHARPERBLOCK  •  512; 

BKGND  -  0;  INK  *  1 ; 

(•  INPUT  COMMANDS  •) 


ON  •  ’T’;  OFF  •  ’F’; 
LINEIN  •  ’L': 

LI ST BLOBS  -  *B* : 

LISTPER IMS  -  ’P* ; 

L ISTRUNLENOTHS  •  ’R’: 
L1STACTIVELINESEGMENTS  - 
DIAGNOSTICS  ■  ’D’ ; 

STOP  •  ’S’; 


(«  THIS  PROGRAM  HANDLES  512  X  512  IMAGES  ONLY 

(•MD*> 

-  NPIXELSPERLINE  /  NBITSPERCHAR  •)  (•MD* ) 

(•MD*> 


’A’; 


•) 


A-6 


TYPE  BINARY  •  BKGND. .  INK 
CNDTYPE  •  SET  OF* 


viwi  lib  -  obi  vr  <•  •  •  Ci  i 

ARRAYS I NT  -  ARRAY  IK.  5]  OF  INTEGER; 


BLOBPTR  -  f BLOBS; 

PTRPERINSECTION  •  t  PER  INSECT  I  ON; 

PERIMPTR  ■  fPERIN; 

BLOBS  -  RECORD 

COLOR:  BINARY; 

COMP:  INTEGER; 

PER  IN  :  PTRPERINSECTION; 
NPER1NPTS:  INTEGER: 

AREA  ,  XNEAN , YMEAN :  REAL; 

KN I N , XMAX ,YM I N.YMAX :  INTEGER; 
PARENT,  NeXT:  BLOBPTR  END; 


PER INSECT  ION  -  RECORD 

LEFT, RIGHT:  PERIMPTR; 

PREV .NEXT :  PTRPERINSECTION  END; 


PERIR  -  RECORD 

LINE, COL :  INTEGER; 
NEXT:  PERINPTR 


END; 


SEGPTR  -  ♦ SEGMENT; 

SEGMENT  •  RECORD 

ST ARTCOL , ENDCOL :  INTEGER; 
BLOB:  BLOBPTR; 

NEXT;  SEGPTR  END; 


PTRR UNLENGTH  •  tRUNLENGTHS; 

RUNLENGTHS  -  RECORD  _ 

ST  ARTCOL, ENDCOL:  INTEGER; 
NEXT :  PTRRUNLENGTH  END ; 


(•  GLOBAL  VARIABLES  •> 

VAR  BLOB . BLOBSDONE :  BLOBPTR; 

RBCYCLEDPTR:  PERIMPTR; 

ARRAYP:  ARRAY  Cl.. II)  OF  PERIMPTR:, 

ACT  I VEL I NE , PREVSEG , CVRRSEG , SEG :  SEGPTR; 

NEWLINE, NEVSEG.LASTTNEVSEG:  PTRRUNLENGTH; 

CMD.ONOFF:  CHAR; 

CMDSET :  CMDTYPE; 

INTERACTIVE:  BOOLEAN; 

L I NENUM , COL , NEWCOMPNUM , NLREAD .NT I MES :  INTEGER; 

TH ACEBLOBS ,  TR  ACEPER I  MS ,  TRACEACT I VEL I NESEGMENTS , TRACERUNLENGTHS :  CHAR 
TRACED UGNOSTICS:  CHAR: 

INPUTIMAGE.DATAFILE.PERIMFILE:  TEXT; 


INPUTNAME, DATA NAME , PER  1 MNAME :  ARRAY!).. 30]  OF  CHAR; 

BIT:  ARRAY! 1 . . NPIXELSPERLINEJ  OF  BINARY; 
I.CHARTOT.LINEGROUP:  INTEGER; 

TYPELINENUMS:  BOOLEAN; 

HANDONLY , BLOBWR ITTEN :  BOOLEAN;  AREATH:  REAL;  IAREATH: 
TRACEPARENT , RECURS :  BOOLEAN; 

NSK1P:  INTEGER; 

(••P*  NEW  PAGE  «> 

PROCEDURE  PAUSE; 

VAR  DUNNYCHAR;  CHAR; 

BEGIN  _ 

VRITELN( 'TYPE  ANY  CHARACTER  TO  CONTINUE’ ) ; 

READLN; 

READ  4  DUNNYCHAR > 

END;  (•  PAUSE  •> 


(•MD*> 


INTEGER; 


PROCEDURE  TYPERUNLENGTHS; 

VAg  V*™'  PTRRUNLENGTH; 

END;  <•  TTPERUVLENGTHS  •) 

PROCEDURE  TYPEACTIVEUHESHniENTS; 

VAR  ALSPtt:  SEGPTR; 

BEGIN 

END;  (•  TYPEACTIVELINESBUmMlS  *> 


A- 7 


(•tf»  NEW  PAGE  • I 

PROCEDURE  WKlTEBLOfHVAR  DEVICE:  TEXT;  VAR  BLOB:  BLOBPTR) 


BEGIN 

END;  <*  VRITE8L0B  •> 


PROCEDURE  TYPEBLOBSfVAR  DEV ICE: TEXT) ; 

VAR  ALSPTR;  SEGPT1; 

BEGIN 

END;  <•  TTPEBLOBS  •) 

NEW  PAGE  •) 

PROCEDURE  PINPOINT (VAR  DIR.Xl ,Y1 ,X2, Y2: 


INTEGER) ; 


(• 


-I 


DIRECTIONS: 

X 

e  t 
2 


-I 

3 


Y  0  4  . 

1  S  6 

e 

7 

BEGIN 

CASE  DIR  OF 

3,2,1:  Y2 

■ 

Yl-l ; 

4,0  ;  Y2 

■ 

Yl; 

5,6,7:  Y2 

■ 

Y 1  ♦  1 

END; 

CASE  DIR  OF 

3,4, S:  X2 

m 

XI- 1 ; 

2,6  :  X2 

n 

XI; 

0,1.7:  X2 

• 

X1U 

END; 

ENP;  <*  DIRPOINT 

•) 

•) 


PROCEDURE  UNPACKS (VAR  XX:  INTEGER;  VAR  DIR:  ARRAYS I NT ) ; 
VAg^,K,XDIV4096,XDIV512»XDIV64tXDIV8:  INTEGER; 

X  :■  -XX; 

XDIV4996  ;•  X  DIV  4096; 

KDIVS12  :•  X  DIV  512; 

XDIV64  X  DIV  64; 

XDIV8  :•  X  DIV  8: 

DIR  1 1 )  :•  X  -  XDtmS: 

DIR(2)  :•  XDIV8  -  XD1V64*8; 

D1RI31  :•  XDIV64  -  XDIV5t2*8: 

D)R(4)  XDIVSI2  -  XDIV4096«8; 

DIR (S)  :•  XDIV4096; 

END;  <•  UNPACKS  •) 

(••P*  NEW  PAGE  •> 

PROCEDURE  UNPACKPERINSIVAR  PTR.UPPTRR, UPPTRL:  PERINPTR) ; 

VAR  DIR:  ARRAYS! NT; 

1,11:  INTEGER; 

BEGIN 

ARRAYPU ) t .LINE  :•  PTRt.LINE; 

ARRAYP1 l ) « .COL  :■  PTRt.COL  -  NPIXELSPERLINE; 

UNPACKS (PTRt .NEXTt .LINE, DIR) ; 

FOR  I  :■  I  TO  S  DO  BEGIN 
1 1  ;  ■  I  ♦  1 : 

DIRPOINT  (DIR(  1 1 ,  ARRAYPf  lit  .COL,  AKJMYPf  fl  *  .LINE, ARRAYPf  I J)  t  .COL, ARRAYPf  I  J )  t  .LINE) ; 
END;  <*  FOR  •) 

UNPACK5(PTRt .NEXTt .COL.DIR) ; 

IF  TRACEDI AGNOSTICS  •  ON  THEN  FOR  I  :•  1  TO  5  DO  WRITELN(DIRI  11:4); 

FOR  I  :*  6  TO  10  DO  BEGIN 
II  :■  1  ♦  1; 

DlRPOINT<DIR( 1-5) , ARRAYPf Ilf .COL, ARRAYPf I ) ♦ .LINE, ARRAYPU 11  ♦ ,COL,ARRAYPI 1 11 1 .LINE) 
END;  (•  FOR  •) 


UPPTKR  ARRAYPU); 

UPPTRL  ARRAYPU)); 

END;  <•  UNPACKPERINS  •) 

(••P*  NEW  PAGE  •) 

FUNCTION  CONNECTED (VAR  POINT! ,P0INT2:  PER1NPTR):  BOOLEAN; 


A-8 


BEGIN 

IF  UBS<POINTIt  .LINR-P0lNT2t  .LINE)  <-  l)  AND 

(ABS(  POINT  1 1  .COL  -P01NT21  .COL  )  <-  1>  THEN  CONNECT®  :*  TRUE 
ELSE  CONNECTED  :■  FALSE; 

END; 


PROCEDURE  VR1TEP0!NTS(VAR  DEVICE:  TEXT; 

VAR  RIGHT  .LEFT;  PE3UMPTR; 

VAR  NPER1NPTS;  INTEGER); 

VAR  PTR , PTR2 , UPPTRR , UPPTRL :  PERINPTR; 

BEGIN 

(•  IF  (TRACED I AGNOSTICS  ■  ON)  AND  RECURS  THEN  BEGIN 
WRITELN l  ’RECURS-TRUEM  ; 

PTR2  ;•  RIGHT; 

WHILE  PTR2  <>  LEFT  DO  BEGIN 

WRITELN(PTR2t  .LlNE:4,m2«  ,C0L:4> ; 

PTR2  :•  PTR2t .NEXT;  END; 

WRITELN<PTR2t.LINE:4.PTR2t.C0L;4);  END;  •) 

PTR  RIGHT; 

NPERINPTS  :•  NPERINPTS  ♦  Is 

WITH  PTRt  DO  WRITELN<DEV|CE,LINE:4,  *  \C0L:4>  ; 

WHILE  PTR  <>  LEFT  DO  BEGIN 

IF  PTR t. NEXT t. COL  >  NPIXELSPERLINE  THEN  BEGIN 

IF  RECURS  THEN  BEGIN  WR I TELN < ’ 2ND  RECURSION1);  PAUSE  END; 
UNPACKPER INS (PTRt .NEXT .UPPTRR, UPPTRL) ; 

(•  IF  TRACED  I  AGNOSTICS  »  ON  THEN  WRITELN ( ’VRITEPOINTS  —  AM; 

IF  TRACED I AGNOSTICS  «  ON  THEN 
RECURS  :•  TRUE;  •) 

(•  PTR2  :•  UPPTRR; 

WHILE  PTR2  <>  UPPTRL  DO  BEGIN 

WRlTELN<DEVICE,PTR2t . LlNE:4,PTR2t .COL: 4) ; 

PTR2  :•  PTR2* . NEXT ;  END; 

WR!TELN(DEV|CE,PTR2t .LINE:4,PTR2t .COL:4) ;  •) 

WR ITEPO I NTS (DEVICE, UPPTRR , UPPTRL . NPER INPTS) : 

(•  IF  TRACED 1 AGNOST 1 CS  *  ON  THEN  WRITELN ( ’ WR1TEPOINTS  —  BM  ; 

IF  TRACEDiAGNOSTICS  •  ON  THEN  RECURS  :«  FALSE;  •) 

PTR  :■  PTR». NEXTt. NEXT;  END 
ELSE  BEGIN 

<*  IF  NOT  CONNECTED < PTR , PTRt . NEXT )  THEN  BEGIN 
WRITELN (’POINTS  NOT  CONNECTED’); 

WH I TELN (DEV ICE, ’POINTS  NOT  CONNECTED M ;  PAUSE;  END;  •) 

(•SKIP  UNNECESSARY  PERINETFR  POINTS  •) 

IF  (PTRt. NEXT  <>  LEFT)  AND  CONNECTED (PTR, PTRt .NEXTt .NEXT) 

THEN  PTR  :•  PTRt. NEXT; 

PTR  PTRt. NEXT; 

NPERINPTS  :•  NPERINPTS  ♦  1; 

(•  IF  TRACEDIAGNOSTICS  «  ON  THEN 

WITH  PTRt  DO  WRITELN(LINE;4, *  \COL;4>;  •) 

WITH  PTRt  DO  WRITELN(DEVICE,LJNE:4/  1 ,C0L:4) ; 

END;  (•  ELSE  •) 

END;  <«  WHILE  *) 

END;  (•  VRITEPOINTS  M 

(♦«*♦  NEW  PAGE  •> 

PROCEDURE  VRITEPERINSlVAR  DEVICE:  TEXT;  VAR  BLOB:  BLOBPTR); 

VAR  PSPTR:  PTRPER INSECT  ION; 

B0GIN 

WITH  BLOBt  DO 

IF  PERIN  <>  NIL  THEN  BEGIN 
NPERINPTS  :•  0; 

IF  NOT  HANDONLY  THEN  WRITELN  (DEV  ICE, 1  Ms 

IF  NOT  HANDONLY  THEN  WRITELN (DEV ICE, ’BLOB  CONPONENT  #  1 ,C0NP:4> 
PSPTR  :•  PERIN; 

REPEAT 

IF  NOT  HANDONLY  THEN  WRITELN I DEV ICE, •  M; 

IF  NOT  HANDONLY  THEN  WRITELN (DEV ICE. ’LINE  COL’): 

WITH  PSPTR  t  DO  WR  ITEPO  I  NTS  (DEV  ICE,  RIGHT,  LEFT,  NPER INPTS) ; 
PSPTR  :■  PSPTRt . NEXT 
UNTIL  PSPTR  •  PERIN; 

END  <•  IF  •) 

END;  (•  VRITEPERINS  •> 


PROCEDURE  TYPEPERINSiVAR  DEV  ICE:  TEXT ) ; 

VAR  ALSPTR :  SEGPTR; 

BEGIN 

WR!TVXN(DEVICE, 1  M; 

WRITELN  (DEV  ICE,  ’BLOB  PERINETEXSM  ; 
ALSPTR  :•  ACTIVELINE: 

WHILE  ALSPTR  <>  NIL  DO  BEGIN 

WRITEPERINSIDEV  ICE, ALSPTRt  .BLOB) ; 
ALSPTR  :•  ALSPTRt. NEXT  END; 

END;  <*  TYPEPERINS  •) 


v 


ME W  PAGE  •> 

PROCEDURE  MEW  (VAX  NEWPO  INTER:  PERIHPTX); 

(•  THIS  ROUTINE  IS  NEEDEO  ONLY  IP  THE  COMMAND  'DISPOSE1 
IN  THIS  VEXSiON  OF  PASCAL  *) 

BEGIN 

IP  BECYCLEDPTX  <>  NIL  THEN  BEGIN 
NEWPO  INTER  :•  RECYCLEDPTR; 

RECYCLEDPTR  :•  RECYCLEDPTR!  .NEXT: 

NEVPOINTEKt .NEXT  :•  NIL  END. 

ELSE  NEW (NEWPO INTER) ; 

END;  (*  NEW  •) 


<»NSP*> 

IS  NOT  AVAILABLE 


PROCEDURE  DELETEP ERINS(VAR  RIGHT, LEFT:  PER l MPT Rf ; 

VAR  POINT:  PERUXTR; 

BEGIV 

(•  ADD  THE  LIST  OF  PERIMETER  POINTS  OP  THE  BLOB  TO  THE  LIST  OF 
RECYCLED  PERIMETER  POINTS,  THIS  SECTION  OF  CODE  SHOULD 
BE  USED  IF  THE  COMMAND  'DISPOSE'  IS  NOT  IMPLEMENTED  IN  THIS  VERSION 
OF  PASCAL  •) 

LEFTt  NEXT  :■  RECYCLEDPTR; 

RECYCLEDPTR  :•  RIGHT; 

(•  DISPOSE  OF  PERIMETER  POINTS  AND  PERIMETER  SECTION.  THIS  SECTION  OP  CODE 
CAN  BE  USED  ONLY  IP  THE  COMMAND  'DISPOSE'  IS  IMPLEMENTED.  •) 

(•  WHILE  RIGHT  <>  LEFT  DO  BEGIN 
POINT  :■  RIGHT; 

RIGHT  :•  RIGHT* .NEXT; 

DISPOSE (POINT)  END; 

DISPOSE  (RIGHT)  •> 

END;  (•  DELETEPERIMS  •> 

NEW  PAGE  •)  _  _ 

FUNCTION  DIRECT ION(VAX  Xl,Yl,X2,Y2:  INTEGER);  INTEGER; 


DIRECTIONS: 

X 

-1  0  1 

-1  3  2  1 

0  4.0 

1  5  6  7 


•) 


VAR  DELTAY:  INTEGER; 

BEGIN 

DELTAY  :•  Y2  -  Yl; 

CASE  X2  -  XI  OF 

- 1 :  CASE  DELTAY  OF 
-1:  DIRECTION 
0:  DIRECTION 
1:  DIRECTION 


•  3; 

■ 

•  5 


END; 


0:  CASE  DELTAY  OF 

-1 ;  DIRECTION  :■  2; 

1:  DIRECTION  :•  6  END; 


I :  CASE  DELTAY  OP 
•  1 :  DIRECTION 
0:  DIRECTION 
1:  DIRECTION 


1; 

0; 

7 


END:  (•  CASE  DELTAY  •) 

(•  DIRECTION  •> 


END 


END; 

(•tP*  NEW  PAGE  •) 

PROCEDURE  CON? ACT  POINTS (VAR  PTR.PTRB:  PER I MPT R;  VAR  I:  INTEGER); 
V^PR£VL1NE,PREVC0L, POSITION, DIR:  INTEGER; 

"WITH  PTRt  DO  BEGIN 

IF  I  •  0  THEN  BEGIN 
PTRB* .LINE  :•  0; 

PTRBt  .COL  *•  0; 

PTRBt. NEXT  :«  NIL  DID 
£  IF  I  <  6  THEN  BEGIN 

IP  I  ■  1  THEN  POSITION  t*  1  ELSE  POSITION  f  POSITION  •  S: 

PTRBt.  LINE  PTRBt. LINE  -  POSITION  •  DIRECT  ION  <PREVOOL,PREVLINB,COL, LIME) ; 

END 
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BE0IN 

IF  I  •  6  TWK  POSITION  :*  1  ELSE  POSITION  :•  POSITION  •  8; 

PTMt.COL  PTWt.OOL  -  POSITION  •  DIIECTIONIPREVOOL.PnVLINB.OOL.LINB)  i 
Dll); 

PREVLINE  :•  LINE; 

PREVCOL  :•  COL; 

END;  (•  VITN  •> 

END;  <•  CO./*  ACT  POINTS  •) 

l*8P*  NEW  PMS  •> 

PNOCEDURE  PKPERINSIVAN  PEKIN:  PTKPERINSECTION) ; 

TAN  PTNjPTKA.PTKB;  PERINPTR; 

BEGIN 

ptr  rant. right,* 

1  :•  0: 

FZRSTG*OUPOF10  :•  TRUE: 

IF  PTR  <>  PERIM* .LEFT  THEN 
REPEAT  BEGIN 

IF  tPTRt.LINE  <  LIN0nm-l>  AND  (PTRt.OOL  >  0) 

_ _ _ _  AND  < PTR*. COL  <*  NPIXELSPERLINE)  THEN  BEGIN 

IF  I  -  0  THEN  BEGIN 
PTRA  :•  PTR; 

NEVWlPTRB)  END: 

CONPACTPO I NTS ( PTR , PTRB , I) ; 

I  :•  I  ♦  1; _ 

IF  I  -  11  THEN  BEGIN 

(•  TO  IDENTIFY  THAT  COMPACTED  RECORDS  FOLLOW  •> 

PTRAt.COL  :•  PTRA t  .COL  ♦  NPIXELSPERUNE; 

<•  INSERT  COUP ACTED  PERIMETER  RECORD  INTO  THE  LIST  •> 

PTRBt  NEXT  :«  PTRt .NEXT ; 

DELETEPERINS(PTRAt .NEXT. PTR) ; 

PTRA t. NEXT  PTRB; 

PTR  PTRB; 

I  :■  0; 

END;  (•  IF  I  ■  H  •) 

END  (•  IF  3  CONDITIONS  THEN  *> 

ELSE  IF  I  >  0  THEN  BEGIN 

<•  AN  INCOMPLETE  COMPACTION  OCCURED  —  RESET  VARIABLES  •> 

I  :•  0; 

FIRSTGROUPOFJ0  TRUE; 

DELETEPER 1 NS ( PTRB . PTRB )  END ; 

PTR  PTRt .NEXT; 

END  <•  REPEAT  BEGIN  •) 

UNTIL  PTR  •  PERlRt.LEFT; 

END;  <*  PAPER IMS  •> 


PROCEDURE  PACKPER I MS ; 

VAR  NPERIMPTS:  INTEGER; 

BEGIN 

SEO  :•  ACTIVELINE; 

WHILE  SEG  O  NIL  DO  BEGIN 
WITH  SEGt.BLOBt  DO 

IF  PERIM  <>  NIL  THEN  BEGIN 
PKPERIMSIPERIM) ; 

PEKIN  :•  PERIM* .NEXT  END; 
SEG  :•  SEGt.NEXT  END; 

END;  (•  PACKPER INS  •> 

(••P*  NEW  PAGE  •) 

FUNCTION  END0F8L0CK: BOOLEAN; 

BEGIN 

IF  CHAJtTOT  •  NCHARPERBLOCK  THEN  BEGIN 
CHARTOT  :■  0; 

ENDOFBLOCX  TRUE  END 
ELSE  ENDOFBLOCX  :•  FALSE; 

CHARTOT  CHARTOT  *  1; 

DID;  <•  ENDOFBLOCX  •) 


PROCEDURE  GETLINE; 

(•  OET  A  LINE  OF  BINARY  IMAGE  DATA  FROM  AN  INPUT  FILE.  •> 

VAR  t NPUTBYTE * BOBL :  CHAR: 

I ,  J* 1NDX.BITS:  INTEGER; 

BEGIN 


m 


<•  READ  IN  ONE  LINE  OP  DATA  FROM  THE  INPUT  PILE  •> 

«•  UNPACK  THE  BINARY  DATA:  IT  IS  ASSUMED  THAT  THE  DATA  IS  INTERPRETED  AS 
A  SERIES  OF  8-BIT  CHARACTERS  •) 

INDX  0; 

FOR  I  :«  1  TO  NCHARSPERUNE  DO  BEGIN 

IF  ENDOFBLOCK  THEN  READ < INPUT l MAGE, EOBL) ;  (•ND#) 

READ! INPUT  IMAGE, 1NPUTBYTE) ; 

BITS  ORDUNPUTBYTE» ; 

FOR  J  ;•  1  TO  NB I TSPERCHAR  DO  BEGIN 
INDX  :*  INDX  *  l; 

IF  ODDCBITSj  THEN  BIT! INDX)  :*  INK 

ELSE  BIT (INDX)  BKGND; 

BITS  :■  BITS  DIV  2; 

END, 

END; 

END;  <*  GETLINE  •) 

PROCEDURE  SETLINE; 

(•  SET  UP  STARTING  RUNLENCTH  RECORD  ♦> 

BEGIN 

IF  LINENUN  -  1  THEN  BEGIN 
NEW (NEWLINE) ; 

NEVLINEf .NEXT  ;•  NIL  END; 

NEWLINE#. ST ARTCOL  :•  NININT* 1 ; 

NEWL I NE# . ENDCOL  :•  MAXINT-l; 

NEWS EG  NEWLINE; 

LASTNEWSEG  :•  NEWLINE; 

END;  (*  SETLINE  *> 


PROCEDURE  ADDRUNLENGTH (COLUMN:  INTEGER); 

VAR  NEWESTSBO:  PTRRUNLENGTH ; 

BEGIN 

IF  NEVSEGf .NEXT  <>  NIL 

THEN  NEWESTSEG  NEWSEG#  .NEXT  <•  ADVANCE  NEW  SEGMENT  POINTER  IF 

A  SEGMENT  ALREADY  EXISTS  •> 

ELSE  BEGIN 

NEW  (NEWESTSEG);  (•  CREATE  A  NEW  SEGMENT  IF  IT  DOESN’T  EXIST  •) 

NEVSEGf  NEXT  :•  NEWESTSEG;  (•  CONNECT  NEWEST  SEGMENT  TO  LIST  *) 

NEWESTSEG t. NEXT  NIL  END; 

NEVSEGf .ENDCOL  :«  COLUMN  -  I; 

NEWESTSEG f .STARTCOL  :«  COLUMN; 

NEWESTSEG# . ENDCOL  :«  MAX I NT- 1; 

NEWS  EG  :»  NEWESTSEG  (•  UPDATE  SEGMENT  POINTER  TO  THE  NEWEST  SEGMENT  •) 

END;  (•  ADDRUNLENGTH  •) 


PROCEDURE  CR EATER l!N LENGTHS ; 

(•  FOR  EACH  TRANSITION,  CREATE  A  NEW  RECORD  FOR  THE  RUNLENGTH  SEGMENT  •> 

VAR  I:  INTEGER; 

BEGIN 

IF  B1T111  -  INK  THEN  ADDRUNLENGTH U ) ; 

FOR  I  2  TO  NPIXELSPERLINE  DO 

IF  BIT! I- 1 1  <>  HIT! 1 1  THEN  ADDRUNLENGTH < l > ; 

IF  BIT  (NPIXELSPERLINE)  -  INK  THEN  ADDRUNLENGTH  (NPIXELSPERLINE*  I ) ; 
LASTNEWSEG  :«  NEWSEG; 

IF  TRACERUNLENCTHS  •  ON  THEN  TYPERUNLENGTHS ; 

END;  (•  CREATERUN LENGTHS  •) 

(••P*  NEW  PAGE  •> 

PROCEDURE  ADDR I GHTPER INPOINT (VAR  PER IN :PTRPER INSECT  ION;  LINE, COL: INTEGER); 

(•  THIS  PROCEDURE  ADDS  ONE  POINT  TO  THE  PERIMETER  LIST.  THE  NEW  POINT  IS 
INSERTED  BEFORE  THE  POINT  DESIGNATED  BY  PER INt. RIGHT  •> 

VAR  NEWPO INTER:  PER I MPT R; 

BOGIN 

NEW  (NEWPO  INTER) ;  (*NSP«> 

NEWPO INTER# .LINE  :•  LINE; 

NEWPO INTER#. COL  :•  COL; 

NEWPOINTER# .NEXT  :•  PERIMt .RIGHT; 

PERIMt .RIGHT  :•  NEWPOINTER: 

END;  (•  ADORIGHTPER1NPOINT  *) 


PROCEDURE  ADDLEFTPER  INPOINT  (VAR  PERIM:PTRPER)RSECTION;  LINE, COL:  INTEGER)  ; 

(•  THIS  PROCEDURE  ADDS  ONE  POINT  TO  THE  PERIMETER  LIST.  THE  NEW  POINT  IS 
INSERTED  AFTER  THE  POINT  DESIGNATED  BY  PERIMt. LEFT  •) 
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C*NSP*> 


VAX  NEWPO INTER:  PERINPTR; 

BEGIN 

NEW  (NEWPO  INTER) ; 

NEWPO  INTER* .LINE  LINE; 

NEWPO INTER* .COL  :•  COL; 

NEWPO INTER* .NEXT  :■  NIL; 

PER  IN*. LEFT*. NEXT  :«  NEWPO INTER; 
PER IN*. LEFT  NEWPO INTER: 

END;  (•  ADDLEFTPERIMPOINT  *> 


<•  CONNECT  PREVIOUS  POINT  AND  NEW  POINT  •) 
(*  RESET  LEFT  POINTER  •) 


(••P*  NEW  PAGE  •) 

PROCEDURE  INSERTSEGMENTtVAR  ST ART COL , ENDCOL : INTEGER) ; 

VAR  SURRBLOB:  BLOBPTR ; 

SUKRCOMP :  INTEGER; 

SURRCOLOR:  BINARY; 

A ,B:  SEGPTR ; 

COL:  INTEGER; 

NEWPERIN:  PTRPER  INSECT  I  ON; 

BEGIN 

IF  TRACED I AGNOSTICS  •  ON  THEN  VR!TELN< ’BEGIN  INSERTION  *>; 
Si)  HR  BLOB  ;•  PREVSEG*  .  BLOB ;  (•  SURROUNDING  BLOB  •) 

SUKRCOMP  SURRBLOB* .COMP ; 

SURRCOLOR  SURRBLOB* .COLOR; 

(•  CREATE  A  NEW  BLOB  FOR  THE  NEW  SEGMENT  •) 


NEW  (BLOB) ; 

WITH  BLOB*  DO  BEGIN 

IF  SURRCOLOR  •  BKGND  THEN  COLOR  INK 

ELSE  COLOR  :•  BKGND; 

NEVCOMPNUM  :•  NEWCOMPNUN  ♦  I; 

COMP  :■  NFWCOMPNUM; 

AREA  :*  0.0; 

KMEAN  ;•  0.0; 

YMFAN  :«  0.0; 

XM  IN  STARTCOL; 

XMAX  ;•  ENDCOL; 

YM IN  :•  LINENUM; 

YMAX  LINENUM; 

PARENT  :*  SURRBLOB; 

NEXT  :«  NIL; 

IF  COLOR  *  BKGND  THEN  PERIM  :•  NIL  (*  NO  PERIMETER  LIST  FOR  BACKGROUND  BLOB  •) 
ELSE  BEGIN 

NEW (PER IN) ; 

PERIM*. PREV  :•  PERIM; 

PERIM*. NEXT  PERIM; 

PEKIN*  .RKJHT  :■  NIL; 

ADDRIGHTPERINPOINTCPERIM^INENUII.STARTOOL)  ;  (•  ADD  FIRST  POINT  •) 

PERIM*. LEFT  :*  PFRIMt .RIGHT: 

FOR  COL  ; *  STARTCOL* 1  TO  ENDCOL  DO 

ADDRIGHTPER1MP01  NT  (PERIM, LINENUM, COL)  (•  ADD  REMAINING  POINTS  •> 

END  (•  ELSE  •> 

END;  (•  WITH  BLOB*  •> 

IF  TRACED! AGNOSTICS  -  ON  THEN  WRITELN< ’DBUG  INSERTSEOMENT  2')  ; 

(•  FORM  A  NEW  PERIMETER  SECTION  FOR  THE  SURROUNDING  BLOB  •) 


WITH  SURRBLOB*  DO 

IF  COLOR  <  >  BKGND  THEN  BEGIN 

(•  CREATE  A  NEW  PERIMETER  SECTION  AND  INSERT  IT  BEFORE  THE 
CURRENT  PERIMETER  SECTION  •) 

NEV<NeVPERIM>i 
NEWPERIN (.NEXT  :•  PERIM; 

NEWPbRIN* .PREV  j-  PERIM*. PREV; 

PER IN*. PREV*. NEXT  :«  NEWPERIN; 

PERIM*. PREV  :*  NEWPERIN; 

<•  MOVE  PART  OF  THE  RIGHT  END  OF  THE  CURRENT  PERIMETER  SECTION  TO  THE 
NEW  PERIMETER  SECTION  •> 


WITH  PERIM*  DO  BEGIN 

NEWPERIM*  .RIGHT  RIGHT;  <•  SET  RIGHT  POINTER  OF  NEW  SECTION  •> 

IF  TRACED  I  AGNOSTICS  •  ON  THEN  VR 1TELN<’ DBUG  INSERTSEOMENT  2A‘>; 

WHILE  (RIGHT*. LINE  <>  LINENUM- 1 )  OR  (RIGHT*. COL  <>  ENDCOL)  _  „  .  k 

DO  BEGIN  IF  TRACED! AGNOSTICS  •  ON  THEN  WRITELN( RIGHT*  .LINE, RIGHT*  .COL, LINENUM- 1 , ENDCOL) ; 

RIGHT  :«  RIGHT*. NEXT  END; 

IF  TRACED  I  AGNOSTICS  ■  ON  THEN  VRITELN  ( ’  DBUG  INSERTSBGNENT  2B’>; 

NEWPERIN*. LEFT  :•  RIGHT;  (•  SET  LEFT  POINTER  OF  NEW  SECTION  •) _  _ ,  ^  _ 

IF  (RIGHT*.  NEXT*.  LINE  ■  LINENUM- 1 )  AND  (RIGHT*  .NEXTt  .COL  >  RIGHT*  .COL)  THEN  RIGHT  :•  RIGHT*. NEXT 
WHILE  (RIGHT*. NEXT*. LINE  »  LINENUM- 1 )  AND  (RIGHT*  .NEXT*  .CUL  >  RIGHT*. COL)  DO  BEGIN 
IF  TRACED! AGNOSTICS  «  ON  THEN  WRITELN (RIGHT*  .LINE.RIGHT*  .COL, LINENUM- 1) ; 

RIGHT  :•  RIGHT*. NEXT;  (•  SET  THE  RIGHT  POINTER  OF  THE  CURRENT 
SECTION  TO  THE  ENDPOINT  OF  THE  LINE 
SEGMENT  ON  THE  PREVIOUS  LINE  THAT  IS 
PART  OF  THE  SURROUNDING  BLOB  •) 


IF  RIGHT  •  NEW PERIN* .LEFT  THEN  BEGIN 
COL  :*  BIGHT*. COL; 

RIGHT  :•  RIGHT*. NEXT; 

ADDR1GHTPERJMP0!NT(PER!N,LINENUN-1  ,OOL> ;  END; 

NEVPER IN*  LEFT*  NEXT  :«  NIL; 

FND;  (•  WITH  PERINt  •) 

PEKIN  N EVP ERIN 

END;  <•  WITH  SURRBLOB*  •> 

(•  INSERT  2  SEGNENTS  IN  THE  ACTIVE  LINE  I  NEED  I ATELY  BEFORE  THE  CURRENT 
SEGNENT  POINTER  *> 

IF  TRACED I AGNOSTICS  ■  ON  THEN  VRITELNI ’DBUG  INSERTSEGNENT  3*>; 
NEVlA);  NEV(B>; 

A* .STARTCOL  STARTCOL; 

Bt.STARTCOL  ENDCOL* 1 ; 

A* .ENDCOL  :•  ENDCOL; 

Bt. ENDCOL  CURRSEO* . STARTCOL  -  i; 


A*  .BLOB 
B*  BLOB 
PREVSEG* 
A*. NEXT 
B* . NEXT 
CURRSEO 


•  BLOB; 

-  SURRBLOB; 
NEXT  A; 

-  B; 

•  CURRSEG ; 

•  A; 


(•  RESET  CURRENT  SEGNENT  POINTER  •) 


IF  TRACEACT 1 VEL I NESEGNENTS  •  ON  THEN  TYPEACT 1 VEL I NESEGMENTS ; 

IF  TRACEPER1NS  -  ON  THEN  BEGIN 

VR1TEPER I NS (TTY .SURRBLOB) ;  VR I TEPER I NS  <  TTY , BLOB )  END; 

END;  <•  INSERTSEGNENT  •) 

<••?*  NEW  PAGE  •) 

PROCEDURE  ADDBLOB (VAR  BLOB. TERNBLOB:  BLOBPTR) ; 

(*  THIS  IS  A  RECURSIVE  PROCEDURE  FOR  ADDING  A  BLOB  TO  THE  END  OP 
THE  LIST  OF  COMPLETED  BLOBS  •) 

BEGIN 

IF  BLOB  -  NIL  THEN  BLOB  TERNBLOB 

ELSE  ADDBLOBfBLOB* .NEXT .TERNBLOB) 

END;  <•  ADDBLOB  •> 


PROCEDURE  RECORD BLOB ( VAR  TERNBLOB:  BLOBPTR); 

BEGIN 

IF  (HANDONLY  AND  (TERNBLOB* .AREA  <  AREATH))  OR  BLOBVRITTEN  THEN  ELSE 
UR ITEPER I NS ( PER l NF I LE , TERNBLOB) ; 

IF  HANDONLY  AND  (NOT  BLOBVRITTEN)  AND  (TERNBLOB* . AREA  >-  AREATH)  THEN  BEGIN 
BIXJBVR1TTEN  :■  TRUE; 

WRITELN( 'BLOB  OF  HAND  IS  DONE* ) ;  END; 

WITH  TERNBLOB* .PERIN*  DO  DELETEPERINS(R IGHT .LEFT) ; 

(•  DISPOSE(TERNBLOB* .PERIN) ;  •)  (•NSP«) 

VR I TEBLOB  <  DAT AF I LE , TERNBLOB ) ; 

ADDBLOB (BLOBSDONE, TERNBLOB) ; 

END;  <•  RECORDBLOB  •> 


(••P*  NEW  PAGE  •) 

PROCEDURE  DELETfcSEGNENT ; 

VAR  TERNBLOB.  R I GHTBLOB .  LEFTBLOB  REPLACBLOB ,  ABSORBS  I  JOB :  BLOBPTR; 
TKNMCONP ,  R  IGHT  CON  P ,  LEFTCONP  REPLACCONP .  ABSORBCONP :  I  NT  EGER ; 
RIUITCOLOR:  BINARY; 

NEW A REA:  REAL; 

OLDPERIN:  PTRPER INSECT  ION ; 

PER INLPOINT ,PER INRPOINT. INTERIORPOINT.NEXTPOINT:  PERIRPTR; 

COL , ABSORBRCOL :  INTEGER; 

BEGIN 

IF  TRACED  I  AGNOSTICS  •  ON  THEN  VRITELN( 'BEGIN  DELETING*); 

TERNBLOB  :•  CURRSEG 9 .BLOB;  (•  TERMINATING  BLOB  ♦) 

TERNCONP  :•  TERNBLOB* .COMP ; 

RKrHTBLOB  :•  CURRSEG*  .NEXT*  .BLOB; 

RIGIITCONP  :»  K I  GHTBLOB*  .COMP ; 

RIOHTCOLOR  :•  R I  GHTBLOB*  .COLOR; 

LEFTBLOB  :*  PREVSEG* . BLOB ; 

LEFTCONP  :»  LEFTBLOB*  .COMP; 

IF  RIGHTCONP  -  0 
THEN  BEGIN 

REPLACBLOB  :•  R I GHTBLOB;  ABSORBBLOB  :•  LEPTBLOB; 

REPLACCONP  :•  RIGHTCONP r  ABSORBCONP  :•  LEFTCONP  END 

ELSE  BEGIN 

REPLACBLOB  :»  LEFTBLOB;  ABSORBBLOB  :•  RIGHT  BLOB; 

REPLACCONP  :•  LEFTCONP:  ABSORBCONP  :■  RIGHTCONP  END; 

IF  TRACED I AGNOSTICS  -  ON  THEN  BEGIN 

VRtTELNt  ’REPLACING  COMPONENT  -  * , REPLACCONP . 

•ABSORBED  COMPONENT  • * .ABSORBCONP > : 
VRITELNCRIGHTCOLOR-’. RIOHTCOLOR , '  RIGHT  COUP  •  *  .IIGHTCOHP) 


<•  CONNECT  THE  ENOS  OF  TWO  PERIMETER  SECTIONS  OF  THE  TERMINATING  BLOB  IF  ITS 
COLOR  IS  NOT  BACKGROUND  •> 

WITH  TERMBLOBt  DO 

IF  COLOR  <  >  BAGNO  THEN  BEGIN 

(•  ADD  INTERIOR  POINTS  ON  THE  BOTTOM  LINE  OF  THE  TERMINATING  BLOB  TO  THE 
LIST  OF  PERIMETER  POINTS  •> 

FOR  COL  ;•  CURRSEGt  .STARTCOL*  1  TO  CURRSEGt .ENDCOL- l  DO 
ADDLEFTPERIMPOINT (PERIM. LI NENUN-  1  ,COL)  ; 

<•  IF  THERE  IS  MORE  THAN  ONE  PERIMETER  POINT  IN  THE  BLOB  AND 

IF  THERE  IS  ONLY  ONE  POINT  ON  THE  BOTTOM  LINE  OF  THE  TERMINATING  BLOB, 

IT  IS  AT  BOTH  THE  LEFT  END  CF  THE  CURRENT  PERIMETER  SECTION  AND  THE 
RIGHT  END  OF  TIIE  NEXT  PERIMETER  SECTION;  DELETE  THE  POINT  AT  THE 
RIGHT  END  OF  THE  NEXT  PERIMETER  SECTION  •> 

IF  (PERIMt .LEFT  <>  PERIMi .RIGHT)  AND  (CURRSEGt .STARTCOL  ■  CURRSEGt . ENDCOL)  THEN 
WITH  PERIMt. NEXTt  DO  RIGHT  :«  RIGHT*. NEXT; 

IF  TRACED  I  AGNOSTICS  •  ON  THEN  VR ITELN < ’DISPOSE  CHECK  1  '  > ; 

IF  TRACEPERIMS  •  ON  THEN  WRlTEPERiMS(TTY.TERMBLOB) ; 

(•  CONNECT  THE  LOT  END  OF  THE  CURRENT  PERIMETER  SECTION  WITH  THE  RIGHT 
END  OF  THE  NEXT  PERIMETER  SECTION  •> 

WITH  PERIMt  DO  BEGIN 

LEFT  t .  NEXT  :•  NEXTt. RIGHT; 

NEXTt. RIGHT  RIGHT;  (•  RESET  THE  RIGHT  END  POINTER  OF  THE 

NEXT  PERIMETER  SECTION  •)  _ 

PREVt .NEXT  ;•  NEXT;  (•  CONNECT  THE  PREVIOUS  SECTION  •> 

NEXTt. PREV  :•  PREV  END;  <•  WITH  THE  NEXT  SECTION  •> 

OLDPER1M  :•  PERIN; 

PERIM  :•  PERIMt. NEXT;  (•  ADVANCE  THE  PERIMETER  SECTION  POINTER  •> 

(•  IF  THERE  IS  MORE  THAN  ONE  PERIMETER  SECTION,  THEN  DISPOSE  OF  THE  OLD  ONE; 
OTHERWISE,  ADD  THE  BLOB  TO  THE  LIST  OF  BLOBS  DONE  •) 

IF  PERIM  <>  OLDPERIN 

THEN  <•  DISPOSE (OLDPERIN)  •)  (•NSP*) 

ELSE  RECORDBLOB(TERMBLOB) ; 

END;  <*  IF*  WITH  ♦> 

IF  TRACED  I  AGNOSTICS  -  ON  THEN  WRITELN( ’DISPOSE  CHECK  2*>; 

IF  TRACEPERIMS  -  ON  THEN  WRITEPER/MS(TTY, TERMBLOB) ; 


(•  CONNECT  THE  ENDS  OF  THE  TWO  PERIMETER  SECTIONS  OF  THE  REPLACING  BLOB  IF  ITS 
COLOR  IS  NOT  BACKGROUND  •) 

VITH  REPUCBLOBt  DO  BEGIN 

IF  COLOR  <  >  BKGND  THEN  BEGIN 

(•  CONNECT  LEFT  END  POINT  OF  ABSORBED  SECTION  TO  A  PERIMETER  POINT  IN  THE 
REPLACING  SECTION  •> 

PER 1 NLPO I NT  :«  PERIMt . RIGHT ; 

WHILE  PER l NLPO I NT t  COL  <>  ABSORBBLOBt .PERIMt .LEFT f .OOL- 1  DO 
PERINLPOINT  :■  PERlMLPOINTt . NEXT s 
ABSORBBLOBt. PER IMt.LEFTt. NEXT  :«  PERINLPOINT; 

IF  TRACED  I  AGNOSTICS  -  ON  THEN  _ 

WR ITELN ( ' PER  I MLPO I  NT  «’ .PERlMLPOINTt ,LJNE:4, *  t ’.PERlMLPOINTt .COL:4> ; 

IF  TRACEPERIMS  •  ON  THEN  VRlTEPERINStTTY .REPLACBLOB) ; 

PERINRPOINT  :«  PERIMt .RIGHT; 

ABSOR8RCOL  :«  ABSORBBLOBt . PERIMt , NEXTt ,R IGHTt .COL; 

IF  PERIMt, RIGHTt. COL  >  ABSORBRCOL+  I  THEN  _ 

WHILE  PER IMRPOINTt  .COL  <>  ABSORBRCOL*  I  DO  PERINRPOINT  PERIMRPOINTt  .NEXT 
ELSE  IF  PERIMt. RIGHTt. COL  <  ABSOKBRCOL- 1  THEN 

FOR  COL  :«  ABSORBRCOL- 1  DOWNTO  PERIMt .RIGHTt .COL* 1  DO 

ADDR IGHTPER I MPO I NT ( ABSOR6BLOB f . PER I M  t . NEXT , L INENUM- J , COL) ; 

IF  TRACED I AGNOSTICS  -  ON  THEN  WRITELN ( ’DISPOSE  DBUG  2B); 

IF  TRACEPERIMS  •  ON  THEN  WRITEPERIMS<TTY, REPLACBLOB) ; 

<•  DELETE  POINTS  IN  THE  MIDDLE  OF  THE  LINE  SEGMENT  CONNECTING  THE  TWO 
BLOBS  WHICH  ARE  NOT  NOW  PERIMETER  POINTS  •> 

IF  TRACED I AGNOSTICS  -  ON  THEN  _ 

VR ITELN <  ’PERINRPOINT  ■’  .PERIMRPOINTt  .HNE:4,  V  .PERIMRPOINTt  C0L;4) ; 

IF  ( PER  I NLPO I NT  <>  PERINRPOINT)  AND  (PERINLPOINT  <>  PERIMRPOINTt .NEXT)  THEN  BEGIN 
NEXTPOINT  .«  PERIMRPOINTt .NEXT ; 

WHILE  NEXTPOINT  <>  PER  I MLPO I NT  DO  BEGIN 
INTERIORPOINT  :•  NEXTPOINT; 

NEXTPOINT  :•  NEXTPOINTt . NEXT r 

(•  DISPOSE!  INTERIORPOINT) ;  •)  END  END;  (•NSP«) 

IF  TRACED  I  AGNOSTICS  •  ON  THEN  WRITELN( ’DISPOSE  DBUG  4'); 

IF  TRACEPERIMS  -  ON  THEN  VR  ITEPER I  MS  (TTY  ,  REPLACBLOB) ; 

<•  CONNECT  THE  RIGHT  END  POINTS  OF  THE  ABSORBED  AND  REPLACING  SECTIONS  •> 
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IF  PERIMLPOINT  •  PERINRPOINT  THBf 

ADDR 1GHTFER  INPOINT  (ABSORBBLOBt  .PERIM* . NETT, LINENUN.PEMIMRPOINTt .COL) 

ELSE  BEGIN 

PERINRPOINT* .  NEXT  :*  ABSORBBLOB*  .PERI  M*  .NEXT*  .RIGHT; 

ABSORBBLOBt  .PERIM*  .NEXT*  .RIGHT  :■  PERINt  .RIGHT;  END; 

PER IN*. RIGHT  ABSORBBLOBt .PER IN* .RIGHT: 

IF  TRACED 1 AGNOSTICS  •  ON  THEN  VR1TELN (* DISPOSE  DBUG  4B* > ; 

<•  MERGE  PERIMETER  SECTION  RINGS  OF  THE  ABSORBED  BLOB  AND  THE  REPLACING  BLOB  •> 


OLDPERIM  :•  ABSORBBLOBt .PER IN; 

IF  (OLDPERIM  •  OLDPEK IMt .NEXT)  AND  (OLDPERIM  •  OLDPERIMt .PREV) 

THEN  BEGIN  (•  DISPOSE (OLDPERIM)  *>  END  <#NSP#) 

ELSE 

IF  REPLACBLOB  <>  ABSORBBLOB 
THEN  BEGIN 

(•  INSERT  PERIMETER  SECTIONS  OF  ABSORBED  BLOB  BEFORE  PER IMt  •) 
IF  TRACED! AGNOSTICS  «  ON  THEN  WR ITELN ( ’DISPOSE  DBUG  7); 
OLDPERIM*  PREV t .NEXT  :•  PERIM; 

OLDPERIMt  .NEXTt  .PREV  PF.K IMt  . PREV  ; 

PER IMt. PREVt. NEXT  :«  OLDPERIMt .NEXT; 

PER IMt. PREV  OU5PER IMt . PREV ; 

PERIM  :*  OLDPERIM#. NEXT; 

(•  DISPOSE (OLDPERIM) ;  •)  (•MSP*) 

END 

ELSE  BEGIN 


(•  ADD  PERIMETER  LIST  TO  THE  BLOB  DESCRIPTOR  OF  THE  TERMINATING  BLOB  (MUCH 
HAS  THE  BACKGROUND  COLOR  IF  REPLACING  BLOB  -  ABSORBED  BLOB)  •) 

(•  INSERT  PERIMETER  SECTIONS  OF  ABSORBED  BLOB  AFTER  PERU*  •> 
IF  TRACED I AGNOSTICS  ■  ON  THEN  VRITELN( 'DISPOSE  DBUG  8’); 

PERINt .NEXTt .PREV  :•  OLDPERIMt .PREV ; 

OLDPERIMt .PREVt .NEXT  PERINt . NEXT ; 

PERIM  :*  PER IMt .NEXT ; 

TERMBLOBt. PERIM  OLDPERIM; 

OLDPERIMt. RIGHT  :•  PER IMLPOlNT ; 

OLDPER IN*. PREV  :«  OLDPERIM; 

OLDPERIMt  NEXT  :«  OLDPERIM; 

RECORDBLOB ( TERMBLOB)  END; 

IF  REPLACBLOB  <>  ABSORBBLOB  THEN  AREA  AREA  ♦  ABSORBBLOBt . AREA ; 

END;  (•  IF  COLOR  <>  BKGND  •) 


END;  (•  WITH  REPLACBLOB «  •) 


IF  ABSORBBLOB  <  >  REPLACBLOB  THEN  BEGtN 

(•  NOTE  —  THE  PARENT  IS  NOT  COMPUTED  CORRECTLY  •> 

IF  TRACEPARENT  THEN  VRITELNC ’LINE  ‘ , L I NENUM :  3 ,  ’  ABSORBBLOB- * , 
ABSORBBLOBt  .COMP:3,  ’  REPLACBLOB*  •  .REPLACBLOBt  .OOMP)  ; 


(•  CHANGE  ALL  INSTANCES  OF  THE  ABSORBED  COMPONENT  IN  THE  ACTIVE  LINE  TO 
THE  REPLACING  COMPONENT  •) 

SEG  ACTIVELINE; 

WHILE  SEG  o  NIL  DO  BEGIN 

IF  SEGt.BLOB  -  ABSORBBLOB  THEN  SEGt.BLOB  :•  REPLACBLOB; 

WITH  SEGt.BLOB*  DO  IF  PARENT  •  ABSORBBLOB  THEN  BEGIN; 

PARENT  :*  REPLACBLOB; 

IF  TRACEPARENT  THEN  WR  ITELN  ( ’REPLACING  PARENT  IN  SEGMENT  FOR  BLOB  f\C0MP:4>;  END; 
SEG  :*  SEGt .NEXT  END; 

(•  CHANGE  ALL  INSTANCES  OF  THE  ABSORBED  COMPONENT  NUMBER  IN  THE 
LIST  OF  BLOBS  DONE  TO  THE  REPLACING  COMPONENT  NUMBER  •> 

END;  <•  IF  ABSORBBLOB  <>  REPLACBLOB  •> 

(•  DISPOSE  (ABSORBBLOB)  •>  (*NSP*> 

IF  TRACED  I  AGNOSTICS  •  ON  THEN  WRITELNi  ’DISPOSE  DBUG  6M; 

<•  DELETE  THE  CURRENT  SEGMENT  AND  THE  SEGMENT  FOLLOWING  THE  CURRENT  ONE  PROM  THE 
ACTIVE  LINE:  SET  THE  CURRENT  SEGMENT  POINTER  TO  THE  FIRST  SEGMENT  AFTER  THE 
DELETED  ONES  •) 

PREV SEG* .NEXT  CURRSEGt . NEXT * . NEXT ; 

(♦  DISPOSE (CURRSEG*  .NEXT) ;  •>  (♦MSP*) 

(•  DISPOSE (CURRSEG) ;  •)  (*NSP*> 

CURRSEG  :*  PREVSEGt . NEXT ; 

IF  TR ACEACT I VEL I NESEGNENTS  -  ON  THEN  TYPEACTIVELINESBGMENTS; 

IF  TRACEBLOBS  •  ON  THEN  TYPEBLOBS(TTY) ; 

IF  TRACEFERINS  •  ON  THEN  TYPEPERIMS(TTY) 

END;  (•  DELETESEGMENT  •) 
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(••P*  HE*  PAGE  •> 

PROCEDURE  UPDATEPER 1 N ( VAR  PER IN : PTRPER INSECT  ION;  CURRST  ARTCOL , CURRENDCOL , 
NEVSTARTCOL, NEVENDCOL:  INTEGER)  ; 

<•  FOR  CASE  3  ONLY  *> 

VAR  COL:  INTEGER; 

BEGIN 

<•  LEFT  SIDE  OF  PERIMETER  SECTION  •> 

IF  CVRKSTARTCOL* 1  <  NEVSTARTCOL  THEN 

FOR  COL  :•  CURRST ARTCOL* I  TO  NEVSTARTCOL- 1  DO 
ADDLEFTPER INPOINT (PER  IN, LINENUN- I ,COU ; 

IF  CURRST ARTCOL-  1  >  NEVSTARTCOL  THEN 

FOR  COL  CURRST ARTCOL- l  DOWNTO  NEVSTARTCOL* 1  DO 

ADDIJEFTPER I NPO I  NT  <  PER  1 N ,  L I NENUN ,  COL  > ; 

ADDLEFTPERINPOINTiPERIN.LINENUM.NEVSTARTCOL) ; 

PER  IN  •  PER  IN* .NEXT ;  (•  ADVANCE  PERI NET ER  SECTION  POINTER  •> 

<•  RIGHT  SIDE  •) 

IF  CURRENDCOL  >  NEVENDCOL  ♦  1  THEN 

FOR  COL  : ■  CURRENDCOL- 1  DOVNTO  NEVENDCOL* l  DO 
ADDKIGIfTPKRINPOINT  (PER  IN,  LI  NENUN- 1  ,COL)  ; 

IF  CURRFNDOOL  <  NEVENDCOL  -  I  THEN 

FOR  COL  : •  CURRENDCOL* I  TO  NEVENDCOL- I  DO 
ADORIGHTPERINPOINT<PERlM, LINENUN, COL) ; 

ADDRI OHTPER INPOINT (PER IN , LI NENUN, NEVENDCOL) ; 

END,  (•  UPDATEPERIN  •> 


(•BP*  NEV  PAGE  •) 

PROCEDURE  CAS1X)PR0CESSING(VAR  CURRSEG:SEGPTR ;  STARTCOL, ENDCOL:  INTEGER) ; 

VAR  NEVAREA,XL1NEAVE:  REAL; 

BEGIN 

IF  TRACED I AGNOSTICS  -  ON  THEN  WRITELN ( 'BEGIN  CASE  3  PROCESSING* > ; 

WITH  CURRSEGt. BLOBt  DO 

IF  CONP  <>  0  THEN  AREA  ;•  AREA  ♦  ENDCOL  -  STARTCOL  ♦  1; 

IF  (CURRSEGt.  BLOB  t.  COLOR  (>  BXGND)  AND  (CIJRRSEG* .  BLOB* .  TWIN  <>  LI  NENUN)  THEN 

UPDATEPER IN (CURRSEGt  .BLOBt  .PER  IN, CURRSEGt .  STARTCOL ,  CURRSEGt .  ENDCOL ,  STARTCOL ,  ENDCOL) 

IF  TRACEBLOBS  «  ON  THEN  VR I TEBLOB (TTY, CURRSEGt .BLOB) ; 

IF  TRACEPERINS  ■  ON  THEN  VR I TEPER I NS (TTY, CURRSEGt .BLOB) 

END;  (•  CASE3PROCESS I NG  *> 

(••P*  NEV  PAGE  •) 

PROCEDURE  PROCESSLINE; 

VAR  ONEZERO:  INTEGER; 

BEGIN 

IF  LINENUN  •  1  THEN  BEGIN 

NEV (BLOB) ;  (•  CREATE  A  BLOB  DESCRIPTOR  FOR  THE  BACKGROUND  •) 

WITH  BLOBt  DO  BEGIN 
COLOR  :•  BKGND ; 

CONP  :•  0; 

PEKIN  :•  NIL; 

AREA  :*  0.0; 

XNEAN  :•  0.0; 

YMF.AN  0.0; 

XNIN  :•  0; 

XNAX  :•  0; 

YMIN  :•  0; 

YNAX  :•  0; 

PARENT  :•  NIL; 

NEXT  :•  NIL  END; 

<•  INITIALIZE  THE  ACTIVE  LINE  •) 

NEV  (ACTIVELINE) ;  <«  THE  FIRST  SEGMENT  CORRESPONDS  TO  THE  BACKGROUND 

WHICH  IS  TO  THE  LEFT  OF  THE  FIRST  COLUMN  •> 

ACT  I  VELINEt  .STARTCOL  :■  NININTM; 

ACTlVELINEt. ENDCOL  :•  NAXINT-2; 

ACT  I  VELINEt.  BLOB  :«  BLOB:  _  _  _ ,  ^ 

NEV(CURRSEG);  <•  THlS  SEGMENT  DESCRIPTOR  CORRESPONDS  TO  THE  BACKGROUND 
WHICH  IS  TO  THE  RIGHT  OF  THE  LAST  COLUMN  •) 

ACTlVELINEt. NEXT  CURRSEG; 

CURRSEGt. STARTCOL  :•  MAX I NT- It 
CURRSEGt. ENDCOL  :•  NAX1NT-I; 

CURRSEGt. BLOB  BLOB; 

CURRSEGt . NEXT  NIL 
END;  <•  IF  LINENUN  •  I  •) 
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IP  L1NENUM  >  NSKIP  THEN  BEGIN 


CURRSEG  :•  ACTIVELINE;  (•  INITIALIZE  THE  CURRENT  SEGMENT  POINTER  •> 
NEVSEG  NEWLINE;  (•  INITIALIZE  THE  NEW  SEGMENT  POINTER  •) 


REPEAT 

WITH  NEVSEGt  DO  BEGIN 
IF  TRACEDI AGNOSTICS 


ON  THEN  BEGIN 


VI1TELN<  ‘CURRSEGt  .STARTCOL  -  ‘ .CURRSEGt  .STARTCOL,  *  CURRSEGt  .ENDOOL  - 
WR ITELN ( * NEVSEGt . ST ARTCOL  -  * .STARTCOL, ’  NEVSEGt .ENDOOL  -  * , ENDCOL) 

IF  CURRSEGt. BLOBt. COLOR  -  INK  THEN  ONEZERO  :•  1 

ELSE  ONEZERO  :•  0: _ 

WHILE  STARTCOL  >  CURRSEGt . ENDOOL « ONEZERO  DO  DELETESBGMENT :  (•  CASE  1  •> 

IF  TRACED  I  AGNOSTICS  -  ON  THEM  WR  ITELN  <  OBVG  PROCESSLINE  I  * ) ; 

IF  ENDCOL  <  CURRSEGt  .STARTCOL -ONEZERO  THEN  I NSERTSEGNENT ( STARTCOL , ENDCOL ) ; 
CASE3P  ROCE5S  I NG  ( Cl/RRSEG ,  STARTCOL .  ENDOOL  >  ; 

CURRSEGt .STARTCOL  STARTCOL; 

CURRSEGt . ENDCOL  ENDCOL; 

END;  (•  WITH  NEVSEGt  •) 

PREVSEG  :*  CURRSEG;  <•  ADVANCE  THE  PREVIOUS  SEGMENT  POINTER  •) 

CURRSEG  CURRSEGt . NEXT ;  <•  ADVANCE  THE  CURRENT  SEGMENT  POINTER  •) 

NEVSEG  NEVSEGt .NEXT  <•  ADVANCE  THE  NEW  SEGMENT  POINTER  •) 

UNTIL  NEVSEG  •  LAST NEVSEGt. NEXT; 

IF  TRACEACT I VEL I NESEGMENTS  -  ON  THEN  TVPEACT 1 VEL 1 NESEGMENTS : 

IF  CURRSEG  *  NIL  THEN  BEGIN  VRITELNtTTY,  'ERROR:  CURRSBG-N JL* > {PAUSE  END; 

WHILE  CURRSEGt. NEXT  <>  NIL  DO  DELETESBGNENT ;  (•  PROCESS  ALL  SEGMENTS  REMAINING 

IN  THE  ACTIVE  LINE  EXCEPT  THE 
LAST  SEGMENT  •) 


.CURRSEGt .ENDCOL) ; 
END; 


(•  CASE  2  •) 


END;  <•  IP  LINDOm  >  NSXIP  •> 
END;  (•  PROCESSLINE  •> 


BEGIN 


NEW  PAGE  •) 


<*  RAIN  PROGRAM  •) 


HANDONLT  :•  TRUE;  BLOSVRITTEN  :■  FALSE; 

WR ITELN < ’TYPE  INPUT  IMAGE  NAME:1); 

HEADUi :  READ*  INPl/TN4NE>  ; 

RESET < INPUT  IMAGE, INPUTNANE) ; 

KEVR I TE I DATAF I LE, ’ DUMMY . DAT * U 
WRITELNt ’TYPE  OUTPUT  PERIMETER  FILE  NAME 
READLN;  READ (PER INNAME) ; 

REVR 1TEIPER IMP I LE , PER I MNAME) ; 

WR  ITELN  ( 'TYPE  THRESHOLD  FOR  BLOB  AREA 
READLN |  READ (1 AREATH);  AREATH 


I  AREATH; 


nuwvn  ncJiw  mntJiiiw  ,  nnuiin  .  -  uuuuti n i 

VRITELNt •PERIMETER  LIST  COMPACTED  EVERY  X  LINES; 


(•KSP*> 

(•NSP*> 

(•MSP*) 


(•NSP#) 

(•NSF«) 


TYPE  X: ' ) ; 


READLN;  READ < L I NEGROUP ) ; 

TRACEPARENT  :■  FALSE; 

VRITELNt  ’TYPE  0  OF  LINES  TO  BE  SKIPPED  AT  THE  TOP  OF  THE  IMAGE:’); 
READLN;  READ(NSKIP) ; 

RECURS  :■  FALSE; 

WR  ITELN  ( ’TYPE  ”1”  FOR  INTERACTIVE  DEBUGGING’): 

VRITELNt  ’  OR  ”U”  FOR  UNINTERRUPTED  PROCESSING:*); 

READLN;  READ(CND) ; 

IF  CMD  •  *  I  *  THEN  INTERACTIVE  :■  TRUE 


ELSE  INTERACTIVE  :•  FALSE; 


VRITELNt  TYPE  ”T”  TO  TYPE  OUT  LINE  NUMBERS:’); 
READLN:  READ(ONOFF) ; 

IF  ONOFF  •  *T*  THEN  TYPELINENUMS  TRUE; 


LINFNUM  :■  0;  NEVCOMPNUM  :•  0;  BLOBSDONE  :•  NIL;  RECYCLEDPTR  :•  NIL; 
CHARTOT  0; 


NEW  ( ARRAY?  ( 11): 

FOR  I  :•  1  TO  10  DO  BEGIN 
NEWtARRAYPI  I  ♦  1 1  > ; 

ARRAY? (Ilf. NEXT  :«  ARRAYP (1*11;  END; 

ARRAY? U lit. NEXT  NIL; 

IF  NOT  INTERACTIVE  THEN  NLREAD  *•  1; 

TRACEBLOBS  :•  OFF;  TRACEPER1MS  :•  OFF;  TRACERUNLENGTHS  !• 
TRACEACT  I  VEL  1  NESEGMENTS  :»  OFF;  TRACED  UGNOCT  ICS  :•  OFF; 
CMDSET  l’L\  ’B\  *P\  ’A*,  ’R*  *  *D\  ’S’  I; 


OFF; 


REPEAT 

IF  INTERACTIVE  THEN  BEGIN 

WR  ITELN  ( ’TYPE  COMMAND:  ”L”  ,  ”B” .  ’  ’P'  \  ”R”  ,  ”A”  ,  ”1 
READLN:  READ (CMD) ; 

IF  (CMD  -  LINEIN)  THEN  BEGIN 

VRITELNt 'TYPE  NUMBER  OP  LIMES  TO  BE  READ: ’ ) j 
READLN;  READ (NLREAD)  END 
ELSE  NLREAD  :•  I  END 
QJ5E  CMD  :■  LINEIN; 


OR  ”S*”>; 
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/»*  : 


IF  CUD  -  L1NEIN  THEN  FOR  NT  IKES  1  TO  NLREAD  DO  BEGIN 
LINENUN  LINENUN  ♦  I; 

IF  TYPELINENUNS  THEN  WRITELN( ’LINE  \ LINENUN t 3) \ 

GETLINE; 

IF  (LINENUN  -  11  OR  (LINENUN  >  NSKIP)  THEN  BEGIN 

IF  LINENUN  NOD  LINEGKOUP  -  6  THEN  PACK PERI NS; 

SETLINE; 

CNEATERUNLENOTRS; 

PROCESSLINE; 

END} 

END 

ELSE  IF  ICW)  <>  STOP)  AND  (CMD  IN  CNDSET)  THEN  BEGIN 
READ(ONOFF) ; 

CASE  CND  OF 

DIAGNOSTICS:  TRACED UGNOST 1 CS  :•  ONOFF; 

LI ST BLOBS:  TRACEBLOBS  ONOFF; 

L1STPERINS:  TRACEPER IMS  ;■  ONOFF; 

LI ST R UNLENGTHS:  TRACER UNLENGTHS  :•  ONOFF: 
LISTACTIVELINESEGNENTS:  TRACEACT I VEL I NESECNENTS  :•  ONOFF 
END;  (*  CASE  •) 

END;  <•  IF  •) 

UNTIL  (CND  -  STOP)  OR  (LINENUN  ■  NLINES)  OR  BLOBVRITTEN; 

LINENUN  :»  LINENUN  ♦  1; 

SETLINE;  (•  PROCESS  A  LINE  WITH  NO  RUN  LENGTHS  TO  •> 

PROCESSLINE;  <•  ENSURE  CONPLETION  OF  BLOB  PROCESSING  *) 

VRITELNt 'LINEBYLIN  —  ALL  DONE’); 


PROGRAM  SCAN 

C  -  Program  to  rood  enforaetted  seqieatial  dot*  flit*  of  Hoad  data 
C  -  that  vara  geaerated  with  prog  LINEBYLIN  oa  E IKON IX  aoaaa, 

C  ~  aad  attract  position  ia  tpaoa  of  tho  f  lagers  aad  tha  oorresposdlag 
C  -  iatardigital  spaces. 

1  FORMAT  (  N,  IX,  14) 

2  FORMAT  </’  T|pt  ia  filaaaaa  of  iapat  data  fit#  *> 

4  FORMAT  <Q,34A1 ) 

5  FORMAT  (Id) 

6  FORMAT  <’  FINGER  »  '  INCLUDES  INTERDIGITAL  SPACES  -  1  *  LITTLE 

♦  FINGER) 

7  FORMAT  (/IX,  ’FINGER  i\3X, ’ARRAY  POSITION’ .3X, ’XV AL’ .5X’ TV  AL’ , 
♦SX,  ’LENGTH’ ,4X,  ’WIDTH*  ,ZX,  ’RATIO  L/W’  ,3X  ’SLOPE’  ,4X,  ’DiST ) 

8  FORMAT  (ZX.IZ.ieX.IS.SX.FS^.X.FS^.X.FS^, 

♦X,F8.2,X,F8.2,X,F8.2,2X,F8.3> 

9  FORMAT  (/!’  Leagth  of  scat  vtsdov  ( l2"6aa)  ■  *) 

10  FORMAT  (12) 

11  FORMAT  (/IX,’  Laa«th  of  soaa  viadov  •  ’,14,’  (12  -  6  as)*) 

12  FORMAT  (/«’  Aai 1  a  for  Hagers  -  ’) 

13  FORMAT  (P6.0) 

14  FORMAT  (/•’  Aegis  for  thaab  -  ’) 

15  FORMAT  </•’  Finger  eagle  -  ' .P6.8, ’  Thaab  eagle  •  \F6.0) 

16  FORMAT  (/////IX, * Iapat  fila  ; , 14A1 ,2X, ’Oa tpat  fila  \14A1> 

17  FORMAT  (2X, I6,6(2X,F8.2) ) 

19  FORMAT  (/’  At  Thaab  Crotoh:  Handvidth  •  ’ ,F6.2,2N, ’Far ia  • 

♦P13. 2,2X, ’ Araa  •  \FI3.2,2X, ’Ratio  <F»*2/A)  ■  \F6.3> 

20  FORMAT  (/#’  Type  i  to  hava  data  priatad  oat  at  aad  of  tost  *> 

21  FORMAT  (/•’  Type  1  to  aark  aaaaara  poiats  oa  pala  plot  *) 

22  FORMAT  (7(X,!5>> 

23  FORMAT  </’  Entire  Head:  Faria  •  \FI3.2,2X, ’Araa  •  ’,F13.2,2X, 

♦  ’Ratio  if  2/A)  -  '  ,P6 .3) 

24  FORMAT  <’  Type  1  if  yoa  vast  to  eheage  paraaotars  *) 

25  FORMAT  (3I3.X,7FI0.3 ,/. < 10X,7F10.3>> 

26  FORMAT  (’  Typo  aimer ic  codas  for  SUBJECT  ID,  SESSION,  TRIAL  ') 

27  FORMAT  (14,14,14) 

29  FORMAT  (/’  Typa  ia  filaaaaa  for  ovtpit  data  filo  ’) 

29  FORMAT  </'  Typa  DIST  aad  R  critarioa  for  oarvatara  aaalM 

30  FORMAT  <F8.2,F8.3> 

31  FORMAT  </’  Typa  aaabar  (odd)  of  poiats  to  skippad  at  start') 

C - - - - - 

DIMENSION  END(2),NID(2),HFJU)<2),TEMP<2) 

DIMENSION  FNGT1PU8,2)AFNGLEN(  19)  ,FNGVID(!8)  .VALUE (30) 

DIMENSION  SLOPE ( 10) ,ZYER0< 10) ,LINE(30> 

DIMENSION  TANGLE < 100) ,DL£NG( 100) ,CURVE< 10) 

INTEGER  1DATAC6000) , FINGER, FNGNUNC 18) , IWIOTH( 18) f ID IN (6600) 
LOGICAL* I  IF!LE(34>, IFIL£2<34) 


tDafaalt  critarioa  sag  to  shift  for  fiagtips 
tDafaalt  oritarioa  aagla  shift  for  thaab 
!Defail t  dist  batvaaa  poiats 
(Starting  dist  f roa  fiagtip  for  oarva 
tCritariof  fit  for  oarvetire  fa. 

(Skip  thaaa  first  poiats 


ANGLE]  •  90.  tDafaalt  oritarioa 

AMGLE2  •  90.  tDafaalt  oritarioa 

IDIS  -  18  tDafaalt  dist  batv 

DISTC  -  40.0  (Starling  dial  fro 

CRIT  -  0.96  (Criterion  fit  for 

ISTART  -  201  (Skip  thaaa  first 
l FLAG  -  0 
IF  -  0 
TYPE  2 

ACCEPT  4,  IS) , ( IFILEC J> ,  J-1,ISI> 
TYPE  28 

ACCEPT  4,  IS2, ( 1FILE2C J> ,  J-I.IS2) 


TYPE  26 

ACCEPT  27,  IAUTH, ISESS, ITV1AL 
TYPE  26 

ACCEPT  5,  IFLAG 
TYPE  21 

ACCEPT  5,  HPLOTF 

TYPE  24 

ACCEPT  5,  IF 

IF  (IF  .BQ.  0)  00  TO  49 

TYPE  9 

ACCEPT  5,  IDIS 
IDIS  -  IDIS/2 
TYPE  12 

ACCEPT  13,  ANGLE  1 
TYPE  14 

ACCEPT  13,  ANGLE2 
TYPE  29 

ACCEPT  30,  DISTC.CRIT 
TYPE  31 

ACCEPT  S.  I ST ART 

49  CALL  ASS I ON ( 1 , XFILE, IS1) 

DO  SO  X* 1  32000  2 

99  READ  <  I ! I ,ENP*51 tERR«5l)  IDIN(K) ,  IDliUK+l) 

STOP  *  ?  NO  EOF  7  * 

St  IF  <!DSN<K>  .HE.  0)  OO  TO  52 

K  -  K  -  2  tDieregard  OTtrru  poiats 

00  TO  St 
52  NPTS  •  X 

CALL  CUOSEU) 

CALL  SORT <  JDATA, JDIN,NPTS* 1) 

C 

£  -  *■#  !•  elgori tha 

C  -  laitialixe  coeaters  -  find  till*  at  eove  throvgh  data,  eoastaatly 
C  -  apdatiag  the  average  eagle.  nova  through  data  with  a  viadov 
C  -  of  aixa  IDIS,  in  steps  of  IDIS/2.  Poiets  evegd.  with  triaagla  fa. 

C  -  Vhea  fiad  a  larga  shift  ia  the  average  eagle  (ANGLE  for  fiagers, 

C  -  ANGLE2  for  thaab)  store  toeatioa  of  tip,  reiaitialixe 
C  -  average  engte.  Store  the  looatioa  of  the  tip  ia  the 
C  -  data  array,  as  well  as  X  aad  Y  poet t toes  ia  apaoe. 

C  -  Calealate  froa  these  both  leagth  aad  width  of  fiagers 
END<2>  -  AVE( 1 DATA, 1 START U> 

ENDU)  -  AVE<  I  DATA,  1ST  ART) 

NIDC2)  •  AVE<  (DATA,  IDIS*  (START*  I  > 

MIDU)  •  AVE< (DATA, IDIS* I ST ART) 

HEAD(2)  -  AVE< IDATA.2* IDIS* 1ST ART* I I 
HEAD(l)  •  AVEUDaTA,2MDIS*ISTART> 

ANGLE  -  ANGLE  1 
FINGER  *  1 
ICNT  -  0 
N  -  1 

ALPHA  •  F ALPHA (END, HEAD) 

SUKALP  -  ALPHA 
AVEALP  «  ALPHA 

DO  200  I«3#IDIS*  1ST  ART, NPTS,  IDIS 
END(l)  •  MIDI!) 

END<2>  -  NID<2) 

NIDU)  -  HEAD ( 1 ) 

M1D<2>  •  HEAD<2> 

HEAD (2)  *  AVE( IDATA, |e|> 

HEADU)  •  AVE UDATA,!> 

IF  (I  .LT.  N-l>  OO  TO  200 
ALPHA  >  FALPHA (END, HEAD) 

IF  (ABS< ALPHA- AVEALP)  .LE.  100.)  00  TO  70 
IF  (AVEALP  .CTT.  0.6)  ALPHA  -  ALPHA  ♦  360. 

IF  (AVEALP  .LT.  0.0)  ALPHA  »  ALPHA  -  360. 

DALFHA  -  ADS (ALPHA  -  AVEALP) 

IF  (FINGER  .LE.  4)  PRINT  17, 1  ,END(2> , ENDU)  ,HEAD<2> , HEADU ) , 
♦ALPHA, AVEALP 

IF  (ICNT  .LT.  6)  GO  TO  99  tNove  at  least  10ea  before  test 

IF  (D ALPHA  .OT.  ANGLE)  00  TO  100 
N  •  N  ♦  1 
ICNT  •  ICNT  ♦  1 
SDNALP  -  SIDIALP  ♦  ALPHA 
AVEALP  -  SUHALP/N 
00  TO  200 


70 

C 

C 


99 


••  Foasd  a  fiager 


& 
C  «" 
C  - 


H  •  I 
ICNT  •  0 
SUNALP  •  ALPHA 
AVEALP  •  ALPHA 
IF  (FINGER  “ 
IF  (FINOR 


•)  ANGLE  •  ANGUS 

__  _  7)  00  TQ  165 

Par  fiagere  fiae  taae  tha  looatioa  of  the  tip 

If  drevisg  a  iiae  throsgh  the  eeater  of  tho  rigger  aad  flsd  the  poiat 
r  iteraeotioa  with  the  tip. 

L  -  IDGLTS<  IDATA,  1— IDIS— )  ,52.  ,-2> 

N  •  IDELTSt IDATA, I -IDIS- l ,52. ,2) 

J  ■  IDELTSt  IDATA,  I -IDIS- 1 ,72.  ,-2> 

X  •  IDELTSt IDATA, l-IDIS-I ,72. ,2) 


A-20 


156 


160 


165 


190 

200 

201 


(Ylh  -  Yrt>/<xm  -  xjk) 
YJK  -  SLOPE<FINGER)*XJK 


LINE(F!NGKR*4-3)-J-l 
LI HE (FINGER* 4-2) *L- 1 
LIHE(FINGER*4-1 )  •H- 1 
L1NE(F1HGER*4)  «K-1 
XAVE1  •  AVE< I DATA, L) 

XAVE2  -  AVE< I DATA, H) 

YAVE1  •  AVEUDATA.L-l) 

YAVE2  •  AVE( IDATA, N- 1 > 

XLH  •  (XAVEI  ♦  XAVE2>/2 
YLH  -  (YAVEl  ♦  YAVE2J/2 
XAVEI  »  AVE( I DATA, J) 

XAVE2  -  AVEUDATA.X) 

YAVEl  •  AVE<tDATA.J-U 
YAVE2  •  AVEUDATA.X- 1) 

XJK  -  (XAVEI  ♦  XAVE21/2 
YJX  -  (YAVEl  ♦  YAVE2>/2 
SLOPE (FINGER) 

Y2ER0 (FINGER) 

D  -  2. 

LI  -  L 
LI  -  LI  ♦  2 

APRIHE  *  FLOAT UDATA(Ll-l))  -  SLOPE(FIHGO)  •  FLOAT (1DATA(L1 ) ) 
IF  (ABStYZERO (FINGER)  -  APR I HE)  .LT.  D>  GO  TO  160 
IF  (LI  .LT.  N)  GO  TO  150 
LI  -  L 
D  •  2*D 
GO  TO  150 
TYPE  13.  D 
FWmVn(FlNOER)  • 

FNGTIP<FINGER,2> 

FNGT  IP (FINGER, 1) 

GO  TO  190 
FNGNUH(FJNGER)  • 

FNGTIP(FINGER,2) 

FHGT  IP (FINGER, I ) 

FINGER  -  FINGER 
IF  (FINGER  .EQ.  10)  GO  TO  201 
CONTINUE 
1DIS  -  IDIS*2 


LI- 1 

-  IDATA(Ll) 

■  IDATA(Ll-l) 

I-IDIS 
•  NID<2) 

«  HIDU) 
l 


fFoaad  all  fiafar* 


C  -  Caloalata  Jaasth  of  filler 
DO  S10  1*1 ,9,2 
4  •  1 

IF  (1  .GT.  6)  i  •  -1 

C  -  PELTAX  •  FN0TIP(I,2)  -  FNGTIP< !♦ 4,2) 

C  -  DELTAY  *  FNGnPCI.I)  -  FNGT!P<I*J,  1>‘ 

C  -  FNGLENU)  •  SQRT(DELTAX**2  ♦  DELtAY**2) 

FNGLENU)  «  DSTNCE  ( I  DATA ,  FNGNUH  <  I )  + 1 ,  FNGNUH  ( I  ♦  J )  ♦  1  > 

510  CONTINUE 
TYPE  10,1 
C 

C  *  Caloalato  width  of  fiafors  mt  a  dlataaaa  DIST  back  fra*  Up. 

DIST  -  92. 

DO  521  I  *  1,9,2 

L  -  IDIXTSUDATA,FNGNUNU>M,DlST,'-2) 

N  •  IDELTS( IDATA .FNGNUHC DM , DIST, 2) 

C  -  DELTAX  •  FLOAT (IDATA (L>  -  IDATA (H>> 

C  -  DELTAY  ■  FLOAT  ( I  DATA  (L- 1)  -  lDATA(H-D) 

C  -  FNGWIDU)  -  SQRT (DELTAX* *2  ♦  DELTAY* *2) 

FNGVIDU)  •  DSTNCE ( IDATA, L,H> 

DIST  •  123. 

IF  (I  .EQ.  7)  DIST  •  92. 

IF  (I  .EQ.  9)  I  «  10 
WIDTH!  1*2)  -  L-I 
JWIDTHUO)  •  H-l 
521  CONTINUE 
TYPE  10,2 
C 

C  -  Caloalata  ©oreatare  of  fiajarttps  fro*  ♦  to  -  DIST  from  Up. 

C  *  First  find  dalta  an#)*  for  taafaats  aroaad  tha  tips,  thaa  fit 
C  -  tha  aa«la  faactioa  with  a  atratfht  lina  aad  atora  slop#  aad  fit. 
C 
C 


S2S 


$20 


531 

530 


DIST  -  DISTC 
DO  538  I* I, 9,2 

CALL  TANANG(  IDATA, TANGLE, DLENG.DIST.FNGNUHU)*  1 , 6*  12, HP) 
CALL  L!NFIT<DLENG,TANGLE,NP.SL6fE2,R) 

IF  (ABS(R)  GT.  CHIT)  GO  TO  529 

DIST  •  DIST  -  2.0 

IF  (DIST  .LT.  4.0)  GO  TO  529 

00  TO  525 

CURVE! I )  -  SL0PE2 

CURVE! !♦ I >  »  DIST 

DIST  -  DISTC 

TYPE  531 ,  R 

FORHAT  (Fft.3) 

CONTINUE 
TYPE  10,3 


A-21 


I 


C  -  Caloalate  Handwidth,  pariaatar  aad  ir«t  91 1 1  oa  in  f roa  thn  tips  of 
C  -  IKo  liltle  and  iadei  fiagers,  ratpoo t  i  vel j . 

L  -  IDELTSi IDATA.FNGNUM! DM, 185. ,-2> 

N  -  »0ELTStlDATA,FNGNUII(7>M,226.,2) 

C  -  DELTA*  •  FLOAT < I DATA <L>  -  IDATA(M) ) 

C  -  DELTAY  -  FLOAT! IDATA1L-1)  -  IDATA(W-l)) 

C  -  HNDWID  -  SORT !DELTAX#*2  ♦  DELTAY**2) 

HNDtilD  •  DSTNCE! IDATA,L,M> 

CALL  ARPRIN( IDATA,L,N.ARCAtPERIN) 

F1NGAR  •  AREA 
FNOPRN  *  PERIN 
IW1DTH<2>  -  L-l 
(WIDTH! 1 1 >  •  M- l 
TYPE  10,4 
C 

C  -  Fiad  area  aad  pariaatar  —  tarainata  haad  froa  poiat  10aa 
C  -  prosiaal  to  the  tip  of  tha  thaab  to  a  poiat  13ea  prosiaal  to  tha 
C  -  tip  of  the  little  fiager. 

L  •  IDELTS  ! I  DATA ,  FNGNTOH  !  1 )  ♦  1 , 267 . ,  -2> 

N  -  1DELTS (IDATA.FNGNUM! 9)  ♦  1 ,205.  ,2) 

IWI0THU)  •  L- 1 
(WIDTH! 14)  -  N-l 

CALL  ARPRJM<1DATA.L,M, AREA. PERM) 

TYPE  10*5 
C 

C  -  Priat  iafo  aad  write  featara  valaas  to  file.  Faataraa  writtea  la  order: 
C  -  leaath,  width*  rat io!L/W) ,haadwidth,f lager  pariaatar  aad  area, 

C  -  ratio  (fiager  A/D*  haad  pariaatar  aad  area,  ratio  (A/P>. 

RATIO  •  <(FNGPRM/20.S)**2>/(FINGAR/420.) 

RAT  102  -  (< PER IH/20. 5) ••2.)/ UREA/420.  > 

IF  < (FLAG  .NE.  1)  00  TO  551 
C  -  PRINT  22, ! 1V1 DTH! I > ,  I«l,7) 

C-  PRINT  22, (IVJUTH!!),  1-8*14) 

PRINT  16,  ( IF1LE!  I ) ,  I«  1 , 14)  ,  ( IFILE2C  I) ,  IM,14) 

PRINT  15*  ANGLE 1 , AN0LE2 
PRINT  11,  1DIS 
PRINT  6 
PRINT  7 

DO  550  I- 1, FINGER- 1,2 

550  PRINT  8,  <  !♦  1 ) /2 . FNGNUN < 1 )  ,FNGTIP(1 ,2)  ,PWTIP<1 , 1)  , 

♦  PNGLEN (I ) / 20 . 5 , FNOW I D ( I ) /20 . 5 , FNGLEN ( I ) /FNOW I D ( I ) . 

♦CURVE ( I ) f  CURVE 1 1  ♦  I ) 

PRINT  19,  HNDWID/20.5.FNGPRM/20.5,FINGAR/420.  .RATIO 
PRINT  23.  PER  I M/29. 5, AREA/420. .RAT 102 

551  CALL  ASSIGN  tl , IF1LE2, IS2) 

DO  COO  1-1  5 

VALUE (l)  -  FNGLEN(f2-l>/29.5 

VALUE!  I ♦S)  •  FNGWID!  l#2**  1  )/20.5 

VALUE! 1 ♦ 10)  -  FNGLEN! 1*2- 1 l/FROW ID ( 1*2-1 > 

600  CONTINUE 

VALUE! 16)  •  HNDV1D/20. 5 
C  VALUE! 17)  «  FNGPRM/20.5 

C  VALUE! 18)  -  FINGAR/420. 

C  VALUE! 19)  «  RATIO 

VALUE! 17)  -  PER IH/20. 5 
VALUE! 18)  -  AREA/420. 

VALUE! 19)  -  RAT  102 
DO  601  (-1,9,2 

VALUE! !♦ 19)  •  CURVE!!) 

VALUE! I*20>  •  CURVE! !♦ 1) 

601  CONTINUE 

NUNFEA  •  29  _  _ 

WRITE  (1,25)  ( IAUTH,  ISESS,  (TRIAL, ! VALUE!  J) ,  >1,NUNFEA)> 

C 

C  -  Plot  haad  pariaatar  aaabariag  point*  for  fiager  tips  aad  latordlgtal 
C  -  apaoas. 

1Y0FST  -  660 
IXOFST  •  100 

i  -  i 

CALL  !NITT<960) 

CALL  NOVABS!  IXOFST^  IDATA<2) ,  IY0R5T-IDATA<  I) ) 

DO  700  1-3.NPTS.2 

CALL  DRVABS<  lXOFSTalDATAClal) ,  lYOFST- IDATA! I > > 

IF  (NPLOTF  .ME.  1)  00  TO  710 
IF  (1  .NE.  lWItmi(J)>  GO  TO  708 
00  TO  711 

C  -  IF  (I  .NE.  LINE! J) )  GO  TO  700 

710  IF  (I  .NE.  FNGmm!J>)  GO  TO  700 

711  ENCODE! 2±  10,  JJ)  J 
CALL  A0UTST(2,Ji> 

J  ■  J  ♦  1 

700  CONTINUE 

C  -  Flat  liaet  that  iataroapt  fiagar  tips 
GO  TO  1002 
DO  1000  X* I  7 

CALL  NOVAfc<  IXOFST,  I  VOFST.TZEMHK)) 

DO  ieei  I ■ 198.900. 198 

IV  •  IFIX<VZERO<K>  ♦  SLOPE<K)*FLOAT<l)> 

CALL  DRVABSt IXOFST* I , I VOFST* IV> 

1 991  CONTINUE 
1699  CONTINUE 

1992  CALL  FINITT(9,7M) 

STOP 


FUNCTION  F ALPHA (END, HEAD) 

C  -  hiotioi  to  oslonlata  tha  i«|U  of  this  sasaast  of  tho  pariaatsr 


C  -  Aaiitt  jo  fro®  0  to  180  th«a  - 180  to  0 
DIMENSION  ENDtlKHEADU) 

DATA  RAD I AN/57. 296/ ,P 1/3. 1416/ 

DELTAX  •  HEAD<2>  -  END(2> 

DELTAY  •  HEAD(l)  -  END(L) 

IF  (DELTAX  .NE.  0)  00  TO  5 

ALPHA  •  PI/2.  •  RADIAN 

IF  (DELTAY  .LT.  0>  ALPHA  •  ALPHA  •  -1.0 

00  TO  100 

5  ALPHA  -  ATAN2 (DELTAY, DELTAX)  •  RADIAN 

100  FALPHA  •  ALPHA 

RETURN 
END 


FUNCTION  1M0VX< IX, NOVX, I , I DAT A, LARGER, ID I NCT) 

C  -  Posotiofl  to  ®o vo  %  distoooo  NO VX  os  poraaotor  froa  polst  IX  is 
C  -  Airootios  ID1RCT  (♦  for  forward  aovo.  -  for  backward) .  LARGER  Is 
C  -  s  flat  to  toll  whothor  tho  sow  X(Y)  ts  yois*  to  bo  Isrtor  or  sasllor 
C  -  thss  IX.  Retsros  (1)  for  tho  bow  X(Y)  point  is  tho  dots  list. 
DIMENSION  I DATA (6000) 

J  •  I 
L  •  2 

ixprin  -  ix  ♦  novx 

IF  UDIRCT  .LT.  0)  L  •  -2 
IF  (LARGER  .LT.  0)  IXPRIN  •  IX  -  NOVX 
IF  (LARGER  .GT.  0)  00  TO  50 
20  J  -  J  ♦  L 

IF  ( I DATA! J)  .LE.  IXPRIN)  GOTO  100 


00  TO  20 

J  -  J  ♦  L 

IF  (IDATA(i)  .OT. 

00  TO  SO 

I NOVX  -  l 

RETURN 

END 


IXPRIN)  00  TO  100 


FUNCTION  IDELTSC IDATA. L,DIST, J) 

C  -  Fssotios  to  find  tho  loostios  of  tho  X.T  pair  is  tho  array  IDATA  that  is 
C  -  tho  distasco  D1ST  froa  tho  X,Y  pair  poiatod  to  by  L.  Tho  dirootios 
C  -  asd  incr»««t  of  starch  throsjh  tho  array  is  dotoraisod  by  J, 

DIMENSION  IDATA (6000) 

1  FORMAT  (X,  ’ERROR  IN  IDELTS’ > 

2  FORMAT  (X, '  IDELTS:  D  -  \F6.0) 

D  -  2.0 

10  I  -  L 

DELTAS  -  0 

XI  •  FLOAT ( I DAT  A  < l ) ) 

Yl  -  FLOAT ( IDATA* t-l) ) 

20  I  •  I  ♦  J 

X2  •  XI 
Y2  -  Yl 

XI  •  FLOATUDATAU)) 

Yl  -  FLOAT ( I DAT A ( 1 ~ I ) ) 

DELTAX  -  X2  -  XI 
DELTAY  -  Y2  -  Yl 

DELTAS  «  DELTAS  ♦  SORT (DELTAX* *2  ♦  DELTAY**2> 

IF  (ABSCDIST  -  DELTAS)  .LT.  D>  00  TO  100 
C  -  IF  (D  .LT.  24.)  00  TO  25 

C  -  TYPE  1 

C  *  STOP 

25  IF  (I  ,OT.  1)  00  TO  38 

D  -  2*D 
00  TO  16 

30  IF  (1  .LT.  6000)  00  TO  20 
D  •  2*D 
00  TO  10 

100  IDELTS  •  I 

RETURN 
END 


FUNCTION  DSTNCEC IDATA, L.N> 

C  -  Fssotios  to  osloslsts  tho  distsneo  batwooa  two  points,  I DATA (L) 

C  *  sad  1DATA(N)  stint  S  poist  triaafitar  avarsainf. 

DIMENSION  IDAtA<6000) 

HAVEL  •  AVE( IDATA, L) 

RAVEN  -  AVEt IDATA, N) 

YAVEL  -  AVE( IDATA, L-l) 

YAVEN  -  AVE(IDATA.N-I) 

DELTAX  -  XAVEL  -  XAVEN 
DELTAY  -  YAVEL  -  YAVEN 
DSTNCE  •  SQRT (DELTAX* *2  ♦  DELTAY**2) 

RETURN 


oo 


r 

I 

I 


c 

c 


FUNCTION  AVE! I DATA* I) 

Fnaotioa  to  do  »  fivo  point  triii|«ltr  M«ri|«  m4  nlvra  ««ii»  as 
AVE. 

DIMENSION  1 DATA (6000) 

FART  I  •  FLOAT! IDATA! 1-4) >•©.  I 
PART2  •  FLOAT ( I DATA ( 1-2) >*0.2 
FART3  •  FLOAT ! I DATA ( I ) ) *9 . 4 
PAIT4  •  FLOAT ! I DATA < I*2>  >*0.2 
FARTS  -  FLOAT! IDATA! 1*4) >»0.1 
AVE  •  FART  1 +PART2+PART3*PART4+PART5 
RETURN 
END 


FUNCTION  IXYGETOX,  J,  I  DATA) 

C  -  Pan© lion  to  find  t  point  with  at  X(Y>  oovivaWat  to  aaothor  point  (IX) 
C  -  J  ii  tho  starting  point  for  tho  atnroh  through  tho  array. 

C  *  If  J>I0  tho  soarch  will  bo  baokwards  throvgh  tho  array. 

C  -  Rotarna  (I)  for  tho  aov  point. 

DIMENSION  I  DATA  < 6600) 

M  •  2 
L  -  2 

IF  !  J  .Of.  10)  L  •  -2 
I  ■  J 
I  ■  I  ♦  L 

IF  (ABSdDATA(l)  -  IX)  .LT.  M)  00  TO  100 
IF  (I  .QT.  I)  00  TO  30 
M  ■  2*11 
00  TO  10 

IF  (I  .LT.  6000)  OO  TO  29 
M  -  2*M 
00  TO  10 
1XYGET  •  I 
RETURN 
END 


10 

20 


100 


SUBROUTINE  AXFRIM! IDATA, L,R,  AREA  .PER IN) 

C  -  Root  in#  to  oaJcoloto  aroa  and  poriaotor  of  a  portion  of  tho  hand  with 
C  -  booadr its  IDATA(L)  to  IDATA(M) . 

DIMENSION  1DATA(60B0) 

AREA  •  0.0 

FER1N  -  0.0 

X2  •  FLOAT!  (DATA ID) 

Y2  -  FLOAT! 1 DATA! L- I ) ) 

DO  20  l*L*2,M,2 
XI  •  X2 
Y1  •  Y2 

X2  -  FLOAT M DATA! I) ) 

Y2  -  FLOAT ( l DATA ( I - 1 ) ) 

DELTAX  -  K2  -  XI 
DELTAY  •  Y2  -  Y1 

FERIN  •  FER1N  ♦  SORT  (DELTAX* *2  ♦  DELTAY”2> 

AREA  -  AREA  ♦  DELTAX  »  I (Yl*Y2>/2.0) 

20  CONTINUE 

DELTAX  *  FLOAT! 1 DATA <L>  -  I DATA (N) ) 

AREA  -  AREA  ♦  DELTAX  •  ( (FLOAT ( I DATA (M- 1 ) ♦ IDATA  <L- 1 > ) )/2. 0) 

RETURN 

END 


SUBROUTINE  SORT! IDATA.  IDIN.  NETS) 

•tiao  to  aort  data  ao  that  it  alvaya  starta  at  tha  lowor  vriat. 


DIMENSION  IDATA!6000) « IDIN(60B0) 
DO  10  I»NFTS.2,-2 

10  IF! IDIN! I )  .BQ.  1)  00  TO  100 

100  1FUNPTS-I)  .Of.  10)  00  TO  200 

DO  110  I- I, NETS 
1)0  IDATA! I )  -  IDIN! I) 

RETUMN 

200  L  •  0 

DO  210  K-l-i'NFTS 
L  -  L  ♦  I 
I DATA (L)  •  IDIN(K) 

210  CONTINUE 

DO  220  X- 1,1-2 
L  •  L  ♦  1 
IDATA(L)  •  IDtN!K> 

220  CONTINUE 

RETURN 


tStart  with  T 


SUBROUTINE  T AN ANG < 1 DAT A, TANGLE, DLENG .DIST. 

♦JSTART , JSTEP , J RANGE, HP) 

C  -  RoiIim  to  otloilttt  (hi  «i|l«  fiaotioi  troiad  •  s«utio«  of  tho 
C  -  pariwtir.  Tho  section  it  -  to  ♦  DIST  fro*  JSTART.  Tho  loogth  of 
C  “  tho  poriaotor  for  ooch  ton|oot  is  J RANGE,  ond  ioagofttt  oro  ooloolotod 
C  "  for  oooh  soctioo  ia  dots  intervals  JSTEP .  NP  it  tho  ooabor  of  onglot 
C  "  ooloolotod.  Anglot  oro  storod  in  TANGLE  ond  tho  oorrotpondi og 
C  -  loogth  in  DLENG. 

DIMENSION  IDATA(6000) ,TANGLE< 100) , DLENG (108) 

DATA  Pl/3. 1416/, RAD1AN/S7. 296/ 

L  •  IDELTSIIDATA, JSTART, DIST, -2) 

N  •  IDELTSIIDATA, JSTART. DIST, 2) 

NP  -  0 
TLENG  •  0. 

DO  100  !*L,M , JSTEP 
NP  •  HP  ♦  1 

DX  -  AVEUDATA.DJRANGE)  -  AVE<  1  DATA,  I) 

DY  •  AVEUDATA,  I^JRANOE-l)  -  AVE(  1DATA,  1-1) 

C  DX  -  FLOAT  UDATAU  ORANGE)  -  IDATA(l>> 

C  DY  •  FLOAT(IDATA(DJRANGE-l>  -  IDATA(J-1>> 

TLENG  -  fS0irnDX**2  ♦  DY**2>>/2.0  ♦  TLENG 

DLENG (NP)  -  TLENG 

IF  (DX  .NE.  0.)  GO  TO  20 

ANGLE  •  PI/2.  •  RADIAN 

IF  (DY  ,LT.  0.0)  ANGLE  ■  ANGLE  •  -1.0 

00  TO  30 

29  ANGLE  -  ATAN2(DY,DX)  •  RADIAN 

30  IF  (NP  .EX).  1)  ANGLE  1  •  ANGLE 

IF  (ABS( ANGLE- ANGLED  .LE.  180.)  00  TO  40 
IF  ( ANGLE  1  .GT.  0.0)  ANGLE  -  ANGLE  ♦  360. 

IF  (ANGLE!  .LT.  0.0)  ANGLE  -  ANGLE  -  360. 

40  Tangle (Np>  -  angle 

ANGLE 1  -  ANGLE 
100  CONTINUE 

RETURN 
END 


C  - 
C  - 
C  - 
C  - 
C  - 
C  - 
C  - 
C  - 
C  - 
C 


100 

C 

c  - 
c 


SUBROUTINE  UNFIT  (DLENG, TANGLE. NP  ,B,R) 

Rootioo  to  giro  a  toast  sqiaras  fit  to  doto  with  o  iito  T«A*BX. 

X  orroy  of  doto  for  iadopoodoat  rariabio 

Y  array  of  doto  for  depondoot  variable 

NP  naaber  of  poirs  of  doto  poioit 

A  Y  intercept 

SIGMAA  standard  deviation  of  A 

B  slopo 

SI GRAB  stondord  deviation  of  B 

R  litoor  correlation  coefficioot 

DIMENSION  TANGLE ( 108) . DLENG ( 100) 

SUM  -  FLOAT (NP) 

SUNX  -  0. 

SUMY  •  0. 

SUMX2  «  0. 

SUMY2  •  e. 

SUNXY  -  e. 

DO  100  W  ,NP 
XI  -  DLENG ( I ) 

Yl  ■  TANGLE! I ) 
stmx  -  SUNK  ♦  XI 
SUMY  •  SUMY  ♦  YI 
SUNX2  ■  SUMX2  ♦  Xf2 
SUMY2  •  SUMY2  ♦  Y1**2 
SUNXY  •  SUNXY  ♦  XI*YI 
CONTINUE 

Calooloto  eoofficioott  ood  ttoodord  dorioliooa 

DELTA  •  SUM  •  SUMX2-SUMX**2 
A  •  (SUMX2*SUMY-SUMX*SUMXY)/DELTA 
B  •  ( SUMXY*SUN -SUMX*SUMY ) /DELTA 
C  •  NP  -  2 

VARNCE  -  (SUMY2»A*#2*SUM»B#*2#SUMX2-2# 
l ( A*SUNY ♦B#SURXY-A4B#SUMX ) ) /C 

SIGMAA  -  S0RT(VaRNCE4SUNX2/DELTA) 

SIGNAB  -  SQRT(VaRNCE*SUM  /DELTA) 

1  •  (SUN*SUNXY-SUNK*SUNY ) /SORT ( DELTA* (Sim*SUMT2-SUMY442) > 

RETURN 

END 
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I  INTRODUCTION 


This  algorithm  is  taken  from  "Image  Processing  Algorithms  for  In¬ 
dustrial  Vision,"  a  report  by  Gerry  Agin  of  SRI  International.  The 
algorithm  is  referred  to  as  "connectivity  analysis."  It  segments  a 
binary  image  into  "blobs"  (connected  regions)  of  the  same  "color" 

(gray  level).  Only  one  pass  is  made  through  the  image,  and  only  one 
line  is  accessed  at  a  time  (hence  the  name  "LINEBYLINE") .  Features  of 
each  blob  are  computed,  such  as  the  area,  center  of  gravity,  bounding 
rectangle,  and  perimeter  points.  This  document  explains  the  concepts 
used  in  LINEBYLINE  and  provides  an  example  of  processing. 


II  RUN-LENGTH  CODING 


The  first  step  in  processing  is  the  run-length  coding  of  the  binary 
image.  Each  line  of  the  binary  image  is  converted  into  a  series  of  run- 
length  segments.  A  "1"  denotes  the  color  of  "object"  and  "0"  denotes 
the  color  of  "background." 

To  illustrate  how  run-length  coding  is  done,  let  us  consider  the 


following  8x8  binary  image: 


Column  number: 


Row  number :  1 

2 

3 

4 

5 

6 

7 

8 


(Note:  "b"  denotes  an  implicit 


12345678 

bbbbbbbbbb 
b  0  0  0  1  1  1  0  0  b 
b  1  1  0  0  1  1  1  0  b 
bOlOllOllb 
bOllllOOlb 
bOlllOOOlb 
bOlOllOllb 
b  0  1  0  0  1  1  1  0  b 
bOOOOOlOOb 
bbbbbbbbbb 

0"  beyond  the  margins  of  the  image.) 
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The  algorithm  for  run-length  coding  examines  each  pixel  of  the  row 
in  turn  and  records  the  column  number  of  any  element  that  differs  from 
its  predecessor.  A  0  is  assumed  to  precede  the  first  element  and  to 
follow  the  last  one.  For  this  reason,  if  the  last  column  of  some  par¬ 
ticular  row  contains  a  1,  a  transition  will  be  recorded  in  column  9. 

In  general,  there  will  be  twice  as  many  column  numbers  recorded  as 
there  are  segments  of  contiguous  Is  in  the  row. 

The  result  of  run-length  coding  our  example  is  as  follows: 


Row  1 

Col.  4,7 

Row  2 

Col.  1,3, 5, 8 

Row  3 

Col.  2, 3, 4, 6, 7, 9 

Row  4 

Col.  2, 6, 8, 9 

Row  5 

Col.  2, 5, 8, 9 

Row  6 

Col.  2, 3, 4, 6, 7, 9 

Row  7 

Col.  2, 3, 5, 8 

Row  8 

Col.  6,7 

Later  in  this  paper,  we  will  refer  to  the  run-length  segments  that 
make  up  a  line  or  row.  These  are  contiguous  sequences  of  pixels  of 
either  color  (0  or  1).  For  example,  row  1  has  exactly  three  run-length 
segments.  The  first  segment  is  all  Os  and  extends  from  minus  infinity 
up  to  (but  not  including)  column  4.  The  second  segment  in  Is  runs  from 
column  4  to  (but  not  including)  column  7.  The  third  and  last  segment 
consists  of  Os  from  column  7  to  plus  infinity.  Because  every  row  must 
start  and  end  with  0,  there  is  always  an  odd  number  of  run-length  seg¬ 
ments  in  any  line. 

The  starting  and  ending  column  numbers  of  each  run-length  segment 
are  stored  in  a  linked  list  for  later  processing  by  the  connectivity- 
analysis  procedures.  When  the  connectivity  analysis  for  all  the  run- 
length  segments  in  a  line  has  been  completed,  the  run-length  segments 
are  deleted.  The  creation  of  run-length  segments  from  a  binary  image 
line  could  be  easily  performed  in  hardware. 


Ill  INTRODUCTION  TO  CONNECTIVITY  ANALYSIS 


The  purpose  of  connectivity  analysis  is  to  separate  a  binary  image 
into  connected  components.  We  call  these  connected  components  blobs . 

A  hole  is  a  special  case  of  a  blob  entirely  surrounded  by  another  blob 
of  the  contrasting  color.  When  connectivity  analysis  is  applied  to  the 
8x8  example  presented  in  the  previous  section,  the  result  contains 
three  components:  the  background  (component  0),  an  "object"  (component  1), 
and  a  hole  (component  2) .  The  three  components  are  shown  below. 

000  +  +  +00  ...111..  ...+  +  +.. 

+  +00  +  +  +0  11.  .111.  +  +..+  +  +. 

0  +  0  +  +  .  +  +  .1.11.11  .  +  .+  +2+  + 

0  +  +  +  +..  +  .1111.  .1  .+  +  +  +22  + 

0  +  +  +...+  .111.  ..1  .+  +  +222  + 

0  +  0  +  +.++  .1.11.11  .  +  .+  +2+  + 

0  +  00  +  +  +  0  .1..111.  .  +  ..+  +  +. 

00000  +  00  . 1..  . +  .. 

All  pixels  in  the  background  are  assigned  the  same  region  number,  even 
though  in  the  8x8  figure  the  background  appears  to  consist  of  four 
separate  regions.  This  is  because  we  assume  the  image  is  embedded  in 
an  infinite  field  of  0s.  Any  region  of  0s  that  touches  the  margin  of 
the  picture  will  therefore  be  classified  as  background. 

The  algorithm  extracts  and  identifies  connected  components,  com¬ 
putes  useful  feature  information,  and  requires  buffering  of  no  more 
than  one  line  of  data.  This  section  describes  the  general  procedure 
for  connectivity  analysis.  Section  IV  gives  some  details  that 
skipped  in  this  section  for  the  purpose  of  explaining  the  process  more 
easily. 

Connectivity  analysis,  as  implemented  here,  makes  use  of  a  data 
structure,  which  we  call  the  active  line,  to  keep  track  of  the  process¬ 
ing.  The  active  line  consist  of 

•  A  linked  list  of  segment  descriptors.  Each  segment  descriptor 
contain  the  starting  and  ending  column  numbers  of  a  run-length 
segment,  the  component  number  to  which  the  segment  belongs,  and 


a  pointer  to  the  next  segment  descriptor.  A  dummy  segment 
descriptor  gives  the  ending  column  number  of  the  final  run- 
length  segment. 

•  A  pointer  to  the  current  segment  being  processed. 

Before  processing  any  data,  the  active  line  is  initialized  to  con¬ 
tain  a  single  segment  descriptor,  representing  the  background.  Before 
processing  any  row  of  run-length  data,  the  pointer  to  the  current  segment 
pointer  must  be  initialized  to  point  to  the  first  segment.  As  we  start 
to  process  our  example  image,  the  state  of  affairs  is  thus: 

Image  so  far:  bbbbbbbbbb 

Active  line: 

Start  column  -infinity  ^infinity 

Component  number  0  dummy 

Current  segment:  I 

In  the  first  row  of  example  data,  transitions  were  found  in  columns 
4  and  7.  This  implies  that  the  row  consists  of  three  run-length  segments 
Os  from  -infinity  to  (but  not  including)  column  4,  Is  from  column  4  to 
(but  not  including)  column  7,  and  Os  from  column  7  to  +infinity.  After 
processing  the  first  row  we  will  want  the  data  structure  to  represent 
those  three  segments. 

One  procedure  will  examine  each  run-length  segment  in  turn.  Each 
segment  will  take  its  turn  as  the  new  segment,  represented  by  its  start¬ 
ing  and  ending  column  numbers.  The  new  segment  must  be  matched  against 
the  partially  completed  analysis  embodied  in  the  current  line.  In  par¬ 
ticular,  we  must  determine  whether  the  new  segment  overlaps  the  current 
segment,  that  is,  the  segment  of  the  active  line  pointed  to  by  the 
current  segment  pointer.  Three  cases  are  possible: 


Case  1: 


The  two  segments  do  not  overlap  because  the  new  segment  is  to  the 
right  of  the  current  one  (that  is,  the  starting  column  number  of 
the  new  segment  is  larger  than  the  ending  column  number  of  the 
current  segment) . 

Current  segment:  XXXXXXXX 

New  segment :  XXXXXXXX 

Case  2: 

The  two  segments  do  not  overlap  because  the  new  segment  is  to  the 
left  of  the  current  one  (that  is,  the  ending  column  number  of  the 
new  segment  is  smaller  than  the  starting  column  number  of  the 
current  segment). 

Current  segment :  XXXXXXXX 

New  segment:  XXXXXXXX 

Case  3: 

Neither  Case  1  nor  Case  2  obtain,  and  the  segments  overlap. 

Current  segment :  XXXXXXXXXXXX 

New  segment :  XXXXXXXXXXX 

Different  actions  must  be  taken  in  each  of  the  three  cases.  We  will 
explain  each  case  as  it  occurs  in  the  course  of  analyzing  our  example 
image. 

The  first  segment  of  the  new  row  goes  from  -infinity  to  4.  The 
current  segment  pointer  points  to  the  first  segment  in  the  active  line, 
which  goes  from  -infinity  to  ^-infinity.  When  these  two  segments  are 
compared,  we  find  that  Case  3  applies:  the  two  segments  overlap. 

Action  on  Case  3: 

Copy  the  start  column  number  from  the  new  segment  to  the  current 
one.  Then  advance  the  current  segment  pointer  to  the  next  segment 
in  the  active  line. 


Copying  the  starting  number  of  the  new  segment  (-infinity)  to  the 
current  segment  results  in  no  change  in  the  data  structure,  because 
-infinity  was  there  to  begin  with.  The  current  segment  pointer  is  ad¬ 
vanced  to  the  next  segment  descriptor  in  the  active-line  list.  After 
processing  the  first  segment  in  the  first  row,  the  active-line  data 
structure  looks  like  this: 
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Image  so  far: 
New  segment: 


bbbbbbbbbb 
b  0  0  0 


Active  line: 

Start  -inf  +inf 

Component  0  dummy 


Current  segment: 


i 


We  now  analyze  the  next  segment.  This  new  segment  runs  from 
columns  4  up  to  (but  not  including)  column  7.  The  current  segment 
pointer  now  points  to  the  dummy  segment  descriptor,  which  starts  at 
♦infinity.  Therefore  Case  2  applies:  The  new  segment  does  not  match 
any  existing  segment  in  the  active  line.  Room  must  be  made  in  the  data 
structure  for  it. 

Action  on  Case  2: 

Insert  two  new  segments  descriptors  (call  them  A  and  B)  in  the 
active  line  before  the  current  segment.  Let  the  current  segment 
pointer  point  to  A.  Copy  the  ending  column  number  of  the  new 
segment  to  the  start  column  number  of  B.  Choose  an  unused  number 
for  a  new  component  and  let  it  be  the  component  number  for  A.  Copy 
the  component  number  from  the  segment  preceding  segment  A,  to  B. 

Now  proceed  as  in  Case  3. 

The  situation  after  processing  the  second  segment  in  the  first  row 
is  thus: 


Image  so  far : 

b  b 

b  b  b  b  b 

b  b  b 

b  0 

0  0 

New  segment: 

111 

Active  line: 

(A) 

(B) 

Start 

-inf 

4 

7 

♦inf 

Component 

0 

1 

0 

dummy 

Current  segment: 

1 

Two  new  segment  descriptors  were  inserted  in  the  active  line.  Descriptor 
A  represents  part  of  a  newly  discovered  blob  in  the  image;  it  is  assigned 
component  number  1.  Descriptor  B  extends  the  background  down  the  right 
side  of  the  new  blob. 


Descriptor  B  takes  its  starting  column  number  from  the  ending 
column  number  of  the  new  segment  (7)  and  its  component  number  from  the 
segment  before  A  (0).  Having  created  a  new  pair  of  segment  descriptors, 
we  can  proceed  to  match  the  new  segment  to  segment  A  by  performing  the 
Case  3  action.  The  starting  column  number  of  the  new  segment  (4)  is 
copied  to  descriptor  A,  and  the  current  segment  pointer  is  advanced  to 
point  to  segment  B. 

The  third  (and  last)  segment  in  the  first  row  runs  from  column  7 
to  +infinity.  This  will  be  matched  against  the  current  segment,  which 
also  runs  from  7  to  +infinity.  This  is  another  Case  3  overlap,  and  the 
data  structure  after  processing  that  new  segment  will  be: 

Image  so  far:  bbbbbbbbbb 

b  0  0  0  1  1  1 


New  segment : 

0  0b 

Active  line: 

Start 

-inf 

4 

7 

+inf 

Component 

0 

1 

0 

dummy 

Current  segment: 

i 

This  concludes  processing  of  the  first  row. 

The  second  row  has  transitions  at  1,  3,  5,  and  8,  or  five  segments 
to  be  processed.  The  current  segment  pointer  is  reset  to  point  to  the 
first  active  segment. 

The  first  segment  of  this  row  runs  from  -infinity  to  1.  When  it 
is  matched  to  the  first  segment  in  the  active  line,  a  Case  3  overlap  is 
discovered.  The  appropriate  action  (extending  the  segment  representing 
the  background)  will  be  performed. 

The  next  segment  runs  from  column  1  to  (but  not  including)  column  3. 
When  it  is  matched  against  the  current  segment  starting  at  column  4, 

Case  2  applies.  The  newly  discovered  blob  is  assigned  component  number  2, 
producing  the  following  situation  after  processing  the  second  segment: 


Image  so  far: 


bbbbbbbbbb 
bOOOlllOOb 
b 

New  segment:  2  2 

Active  line:  (A)  (B) 

Start  -inf  1  3  4  7  +inf 

Component  0  2  0  1  0  dummy 

Current  segment:  ! 


The  last  three  segments  on  this  second  row  all  produce  Case  3  over¬ 
laps  when  they  are  matched  against  the  active  line.  At  the  conclusion 
of  processing  the  second  row,  the  following  data  structure  obtains: 


Image  so  far:  bbbbbbbbbb 

b  0  0  0  1  1  1  0  0  b 
b  2  2  0  0  1  1  1 

New  segment :  0  b 

Active  line: 

Start  -inf  1  3  5  8  +inf 

Component  0  2  0  1  0  dummy 

Current  segment:  I 


Row  3  has  transitions  at  columns  2,  3,  4,  6,  7,  and  9.  Matching 
the  seven  segments  of  this  row  will  result  in  four  Case  3s,  a  Case  2, 
and  two  more  Case  3s.  After  processing  the  third  row  the  following 
will  be  the  state  of  the  data  structure: 


Image  so  far:  bbbbbbbbbb 

bOOOlllOOb 
b  2  2  0  0  1  1  1  0  b 
b  0  2  0  1  1  3  1  1  b 


Active  line: 
Start 

-inf 

2 

3 

4 

6 

7 

9 

+inf 

Component 

0 

2 

0 

1 

3 

1 

0 

dummy 

Current  segment: 
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In  the  fourth  row  there  are  transitions  at  columns  2,  6,  8,  and  9, 
for  a  total  of  five  segments.  The  first  two  of  these  result  in  Case  3 
overlaps.  But  when  we  match  the  third  segment  of  the  fifth  row  (which 
runs  from  column  6  to  8)  with  the  next  segment  of  the  active  line 
(running  from  column  3  to  column  4)  we  discover  that  Case  1  holds.  Here 
is  the  situation  as  we  discover  the  fact: 


Image  so  far:  bbbbbbbbbb 

bOOOlllOOb 
b  2  2  0  0  1  1  1  0  b 
b  0  2<0>1  1(3)1  1  b 
b  0(2  2  2  2) 


New  segment: 

[0  0] 

Active  line: 

(A) 

(B) 

Start  -inf 

2 

3 

4 

6 

7 

9 

+inf 

Component  0 

2 

0 

1 

3 

1 

0 

dummy 

Current  segment : 

I 

The  current  segment  in 

the  active 

line, 

marked  by 

M<  >.» 

above,  is 

matched  by  any  segment 

in  the 

new 

row. 

The 

new  segment. 

marked  by 

"[  ]”  must  eventually  match  the  segment  denoted  by  "(  )".  We  must 
delete  the  unmatched  segment  "<  >"  from  the  active  line.  Furthermore, 
the  previous  segment,  marked  by  "(  )",  forms  a  "bridge"  between  com¬ 
ponents  1  and  2.  We  must  merge  the  two  components. 

Action  on  Case  1: 

Consider  the  segment  in  the  active  line  before  the  current  one 
(call  it  "A")  and  the  segment  after  the  current  one  (call  it 
"B").  If  the  component  numbers  of  A  and  B  are  different,  then 
merge  the  two  components  by  changing  all  instances  of  B's  component 
number  in  the  active  line  to  A’s  component  number.  Delete  the 
current  segment  and  segment  B  from  the  active  line.  Let  the  cur¬ 
rent  segment  pointer  point  to  the  segment  after  B.  Now  go  back 
to  the  beginning  of  the  entire  matching  procedure,  matching  the 
new  segment  against  the  updated  active  line. 

By  this  procedure,  we  merge  components  1  and  2.  The  merging  changes 
the  component  numbers  of  the  fourth  and  sixth  segment  descriptors  in  the 
active  line  from  1  to  2.  The  current  segment  and  segment  B  get  deleted, 
and  the  current  segment  pointer  is  advanced  to  point  to  the  segment 


starting  at  column  6.  After  we  have  taken  these  steps,  the  data  struc- 
ture  looks  as  follows: 


b  b  b  b 
2  0  0  b 

2  2  0  b 

3  2  2  b 

3  3 


Image  so  far: 

New  segment: 

Active  line: 

Start 

Component 

Current  segment : 


b  b  b  b  b  b 
b  0  0  0  2  2 
b  2  2  0  0  2 
b  0  2  0  2  2 
b  0  2  2  2  2 


(A) 

-inf  2  6 

0  2  3 

I 


7  9  +inf 

2  0  dummy 


Now  the  new  segment  may  be  matched  anew  against  the  active  line  and 


current  segment.  This  time  around,  the  match  is  Case  3,  as  are  the  re¬ 
mainder  of  matches  in  this  fourth  row.  After  finishing  this  row,  the 
state  of  the  data  structure  is  as  shown  below: 


Image  so  far:  bbbbbbbbbb 

b  0  0  0  2  2  2  0  0  b 
b  2  2  0  0  2  2  2  0  b 
b*  02022322b 
b  0  2  2  2  2  3  3  2  b 


Active  line: 

Start  -inf  2 

Component :  0  2 

Current  segment: 


6 

3 


8 

2 


9 

0 


+inf 

dummy 

I 


In  the  fifth  row,  there  are  five  segments  with  transitions  at 


columns  2,  3,  8,  and  9.  Only  Case  3  overlaps  occur. 


The  "image  so  far"  is  not  part  of  the  data  structure;  it  only  is  pre¬ 
sented  as  an  aid  in  visualization  of  the  connectivity-analysis  process. 
We  have  changed  all  Is  in  the  image  to  2s  to  indicate,  conceptually,  the 
merging  of  the  two  components. 


The  sixth  row  has  transitions  at  columns  2,  3,  4,  6,  7,  and  9.  The 
third  segment  (from  columns  3  to  4)  will  generate  a  Case  2  match,  the 
others  will  all  be  Case  3.  The  new  component  created  by  the  Case  2 
match  is  assigned  number  4.  Here  is  what  the  data  structure  looks  like 
after  processing  that  row: 


Image  so  far:  bbbbbbbbbb 

b  0  0  0  2  2  2  0  0  b 

b  2  2  0  0  2  2  2  0  b 

b  0  2  0  2  2  3  2  2  b 

b  0  2  2  2  2  3  3  2  b 

b  0  2  2  2  3  3  3  2  b 

b  0  2  4  2  2  3  2  2  b 

Active  line: 

Start  -inf  2  3  4  6  7  9  +inf 

Component  0  2  4  2  3  2  0  dummy 

Current  segment:  I 


In  processing  row  7,  with  transitions  at  2,  3,  5,  and  8  a  Case  1 
situation  occurs.  The  analysis  of  component  3  is  complete:  Its  segment 
descriptor  disappears  from  the  active-line  data  structure.  The  two 
segments  on  either  side  of  component  3  already  have  the  same  component 
numbers,  so  no  merging  is  necessary: 


Image  so  far:  bbbbbbbbbb 

b  0  0  0  2  2  2  0  0  b 

b  2  2  0  0  2  2  2  0  b 

b  0  2  0  2  2  3  2  2  b 

b  0  2  2  2  2  3  3  2  b 

b  0  2  2  2  3  3  3  2  b 

b  0  2  4  2  2  3  2  2  b 

b  0  2  4  4  2  2  2  0  b 

Active  line: 

Start  -inf  2  3  5  8  +inf 

Component  0  2  4  2  0  dummy 

Current  segment:  1 
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In  row  8,  there  are  transitions  at  columns  6  and  7.  Another  Case  1 
situation  occurs  here,  merging  component  4  with  component  0  (the  back¬ 
ground).  Here  is  the  situation  after  we  have  processed  the  eighth  row: 

Image  so  far:  bbbbbbbbbb 

b  0  0  0  2  2  2  0  0  b 

b  2  2  0  0  2  2  2  0  b 

b  0  2  0  2  2  3  2  2  b 

b  0  2  2  2  2  3  3  2  b 

b  0  2  2  2  3  3  3  2  b 

b02022322b 
b  0  2  0  0  2  2  2  0  b 

b  0  0  0  0  0  2  0  0  b 

Active  line: 

Start  -inf  6  7  +inf 

Component  0  2  0  dummy 

Current  segment :  I 

After  the  last  row  of  image  data,  another  row  of  all  Os  will  com¬ 
plete  the  analysis,  removing  component  2  from  the  active  line. 

We  have  engaged  in  a  little  oversimplification  in  order  to  present 
the  basic  ideas  behind  connectivity  analysis  as  plainly  as  possible. 

In  the  next  section  we  will  fix  the  algorithm  to  operate  correctly  in 
all  cases.  But  before  reading  on,  turn  to  the  beginning  of  this  section, 
where  we  presented  the  three  components  we  wish  the  algorithm  to  extract. 
You  may  verify  that  the  results  we  have  obtained  are  what  was  intended 
except  that  different  numbers  were  assigned. 


IV  CONNECTIVITY  ANALYSIS  IN  DETAIL 


The  simplified  algorithm  presented  in  the  previous  section  is  de¬ 
ficient.  It  has  errors,  and  it  doesn’t  compute  any  features  of  the  blobs 
or  any  perimeter  lists.  In  addition,  processing  of  some  borderline 
cases  of  overlapping  segments  was  not  adequately  described.  At  the  end 
of  this  section  the  complete  algorithm  with  its  amendments  and  corrections 
will  be  presented. 
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A.  Handling  Some  Special  Cases 


Let  us  finish  the  processing  of  our  test  image.  To  complete  the 
analysis,  a  last  row  of  Os  running  from  -infinity  to  ^infinity  is  added 
to  the  active-line  data  structure,  leaving  us  with  the  following  situa¬ 
tion: 


Image  so  far: 


New  segment: 

Active  line: 

Start 

Component 

Current  segment: 


b  b  b  b  b  b  b 
b  0  0  0  2  2  2 
b  2  2  0  0  2  2 
b  0  2  0  2  2  3 
b  0  2  2  2  2  3 
b  0  2  2  2  3  3 
b  0  2  0  2  2  3 
b  0  2  0  0  2  2 
b  0  0  0  0  0  2 
b  b  b  b  b  b  b 


-inf  6 

0  2 

I 


b  b  b 
0  0b 
2  0b 

2  2b 

3  2b 
3  2b 
2  2b 
2  0b 
0  0b 
b  b  b 


7  +inf 

0  dummy 


We  have  just  processed  a  line  with  nothing  in  it,  but  the  active 
line  incorrectly  points  to  a  segment  there.  Tt  takes  a  Case  1  match 
between  segments  to  delete  component  2  from  the  data  structure,  but 
there  is  nothing  left  to  process.  What  we  must  do  is  to  check  at  the 
end  of  each  line  whether  the  current  segment  pointer  points  to  the  dummy 
segment  at  the  end.  If  it  does  not,  then  the  Case  1  processing  must  be 
performed  until  the  current  segment  pointer  does  point  to  the  last  seg¬ 
ment. 


A  second  error  concerns  the  merging  of  two  regions  when  a  "bridge" 
connects  them.  The  rule  for  Case  1  states  that  if  the  two  component 
numbers  are  different,  then  the  number  on  the  right  replaces  the  number 
on  the  left  everywhere  else  in  the  active  line.  However,  if  the  com¬ 
ponent  number  on  the  right  happens  to  be  the  background,  the  replacement 
would  renumber  the  background.  We  wish  the  background  to  be  always 
component  number  0  and  must  modify  the  renumbering  rule  accordingly. 
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The  main  reason  for  doing  connectivity  analysis  is  to  obtain  infor¬ 
mation  about  the  components  we  find  thereby.  The  information  we  have 
derived  related  to  color,  location,  area,  perimeter,  and  points  to  other 
components  that  surround,  neighbor,  or  form  holes  in  the  component  we 
are  examining.  The  algorithms  for  deriving  such  information  will  be 
described  later  in  this  report.  Such  information  is  stored  in  a  blob 
descriptor  for  each  component. 

We  must  modify  our  algorithm  description  to  include  the  use  of  blob 
descriptors.  When  a  new  component  is  started  as  a  result  of  Case  2 
processing,  a  new  blob  descriptor  must  be  created.  When  segments  are 
removed  from  the  active  line  in  Case  1  processing,  we  must  keep  track 
of  the  blob  descriptors  that  are  removed.  Several  segment  descriptors 
in  the  active-line  data  structure  can  have  the  same  component  number; 
hence  they  can  point  to  the  same  blob  descriptor.  Therefore,  when  a 
segment  descriptor  is  removed  as  the  result  of  Case  1  processing,  the 
algorithm  must  check  to  see  if  other  segment  descriptors  in  the  active 
line  refer  to  the  same  component.  If  this  is  the  last  reference  to 
that  component,  its  processing  is  finished. 

At  this  point,  the  blob  descriptor  describes  a  connected  component 
whose  connectivity  has  been  completely  determined.  Furthermore,  the 
extraction  of  blob  features  and  perimeter  lists  is  complete.  We  have 
more  than  one  option  as  to  what  to  do  with  the  blob  descriptor.  We  can 
pass  the  address  of  the  completed  blob  descriptor  to  an  application- 
dependent  subroutine  that  will  do  further  processing,  take  some  action 
depending  on  what  was  found,  or  delete  the  blob  descriptor  to  make  its 
memory  space  available  for  other  blobs  or  other  purposes.  Alternatively, 
the  address  of  the  blob  descriptor  can  simply  be  added  to  a  list  of 
isolated  blobs  in  the  image  and  further  processing  will  be  deferred 
until  the  entire  picture  is  analyzed.  In  the  current  implementation  the 
features  and  perimeter  lists  of  the  blob  are  written  to  output  files, 
and  the  blob  descriptor  is  recycled  to  save  memory  space. 


C.  Treatment  of  Diagonally  Adjacent  Pixels 

Two  pixels  belong  to  the  same  component  if  and  only  if  there  exists 
a  path,  along  adjacent  elements  of  the  same  color,  from  one  pixel  to  the 
other.  When  dealing  with  points  on  a  rectangular  grid,  it  is  reasonable 
to  ask  whether  diagonal  points  are  adjacent*  If  diagonal  points  are  to 
be  considered  adjacent,  then  the  components  of  an  image  are  said  to  be 
8-connected ,  since  each  pixel  has  eight  connected  neighbors.  If  diagonal 
points  are  not  considered  adjacent,  then  the  components  are  called  4- 
connected.  Our  connectivity  algorithm,  as  presented  so  far,  finds  the 
4-connected  components  of  an  image.  If  we  were  to  analyze  the  image  of 
a  checkerboard,  then  each  square  of  the  image,  black  or  white,  would  be 
a  separate  component. 

It  would  be  desirable  to  modify  the  connectivity  algorithm  so  it 
has  the  property  that  every  component  (except  the  background)  is  en¬ 
tirely  enclosed  by  one  and  only  one  component,  which  is  of  the  opposite 
color.  Then  every  boundary  between  two  regions  is  a  closed  curve  that 
separates  an  enclosing  region  from  an  enclosed  region.  The  components 
of  a  scene  thus  analyzed  may  be  arranged  hierarchically,  each  component 
having  a  single  superior  (encloser)  and  zero  or  more  inferiors  (holes 
or  enclosed  regions) . 

A  way  of  achieving  this  is  to  specify  that  all  white  cells  (Is) 
are  8-connected  and  that  all  black  cells  (Os)  are  4-connected.  With 
such  a  convention,  the  image  of  a  checkerboard  would  be  analyzed  as  a 
number  of  individual  white  squares,  all  embedded  in  a  single  black 
region  (the  background) .  To  implement  this  requires  modifying  the 
tests  for  overlap  that  determine  Case  1,  2,  or  3.  The  modification  has 
to  be  such  that  one  set  of  tests  applies  when  a  segment  of  Is  is  being 
considered  and  a  different  set  of  tests  applies  for  a  segment  of  Os. 

D.  The  Connectivity  Algorithm 

We  shall  now  consolidate  what  has  gone  before.  The  algorithm 
will  be  presented  in  an  informal  fashion.  Operations  such  as  creating 
and  modifying  lists,  or  obtaining  storage  for  a  descriptor  are  imple¬ 
mentation  details  that  will  not  be  mentioned  here. 
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First,  some  definitions. 


•  Segment  descriptor  is  an  array  or  a  block  of  data  that  contains 
four  items:  starting  and  ending  column  numbers,  a  pointer  to  a 
blob  descriptor,  and  a  pointer  to  the  next  segment  descriptor 
in  the  active  line. 

•  The  active  line  is  a  linked  list  of  segment  descriptors.  The 
current  segment  pointer  points  to  a  segment  descriptor  on  this 
list.  The  active  segment  is  the  segment  descriptor  pointed  to 
by  the  current  segment  pointer. 

•  A  blob  descriptor  is  a  record  of  data  that  contains  at  least  one 
item:  a  color,  which  may  be  either  0  or  1.  Additional  itmes 

in  a  blob  can  be  used  for  feature  analysis,  as  will  be  described 
later.  A  component  number  is  a  number  that  identifies  a  blob. 

To  process  an  image: 

•  Obtain  a  blob  descriptor  to  represent  the  background.  Set  the 
color  word  of  the  background  blob  to  0. 

•  Initialize  the  active  line  to  contain  two  segment  descriptors. 

The  first  segment  descriptor  should  have  a  column  number  smaller 
than  zero  and  a  component  number  pointing  to  the  background 
blob.  The  column  number  of  the  second  segment  descriptor  should 
be  a  large  positive  number,  a  number  greater  than  the  number 

of  columns  in  the  image.  The  component  number  of  the  second 
descriptor  is  irrelevant. 

•  Process  each  row  of  the  image,  as  described  below. 

•  Finish  by  processing  an  extra  row  consisting  of  all  zeros. 

To  process  a  row: 

•  Initialize  the  current  segment  pointer  to  point  to  the  first 
segment  descriptor  in  the  active  line. 

•  Obtain  the  run-length  representation  of  the  row  (see  Section 
II).  The  run-length  data  should  start  with  a  negative  number 
and  end  with  a  large  positive  number. 

•  For  every  pair  of  adjacent  numbers  in  the  run-length  data,  in 
turn,  call  the  segment-processing  operation  (defined  below). 

* 

•  While  the  current  segment  pointer  does  not  point  to  the  last 
segment  in  the  active  line,  perform  the  deletion  operation  (de¬ 
fined  below).  Repeat  this  step  until  the  current  segment  pointer 
does  point  to  the  last  segment. 


In  this  context,  the  word  while  refers  to  a  condition  to  check  and  an 
action  to  perform.  It  means  to  repeat  the  action  zero  or  more  times 
until  the  condition  is  not  true. 
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To  process  a  new  segment ,  given  a  starting  column  number  and  an 
ending  column  number : 

•  While  the  starting  column  number  is  greater  than  the  ending 
column  number  of  the  current  segment  in  the  active  line,  do  the 
deletion  operation  (defined  below).  (This  is  Case  1.) 

•  If  the  ending  column  number  is  less  than  the  starting  column 
number  of  the  current  segment,  perform  the  insertion  operation 
(defined  below),  passing  on  the  starting  and  ending  numbers  of 
the  segment.  (This  is  Case  2.) 

•  Do  feature  extraction  for  Case  3  processing. 

•  Copy  the  starting  and  ending  column  numbers  of  the  new  segment 

to  the  starting  and  ending  column  numbers  of  the  current  segment. 

•  Advance  the  current  segment  pointer  to  point  to  the  next  segment 
descriptor  in  the  active  line. 

To  perform  the  insertion  operation,  given  starting  and  ending 
column  numbers: 

•  Obtain  the  component  number  of  the  segment  descriptor  preceding 
the  current  segment.  Call  that  component  the  surrounding  com¬ 
ponent  . 

•  Obtain  a  new  blob  descriptor.  Call  it  the  new  component.  Set 
the  color  word  of  the  new  component  to  the  opposite  of  the  color 
of  the  surrounding  component. 

•  Obtain  two  new  segment  descriptors  and  insert  them  in  the  active 
line  immediately  before  the  current  segment,  calling  the  first 
segment  "A"  and  the  second  "B.”  Segment  A  receives  the  new  com¬ 
ponent  number  and  the  starting  column  number.  Segment  B  re¬ 
ceives  the  surrounding  component  number  and  the  ending  column 
number . 

To  perform  the  deletion  operation: 


•  Call  the  component  number  of  the  current  segment  the  terminated 
component .  Call  the  component  number  of  the  segment  preceding 
the  current  one  the  left  component  and  the  component  number  of 
the  segment  following  the  current  one  the  right  component.  Call 
the  left  component  the  replacing  component  and  the  right  component 
the  replaced  component.  If  the  right  component  points  to  the 
background,  then  call  it  the  replacing  component  and  the  left 
component  the  replaced  component;  otherwise,  call  the  left 
component  replacing  and  the  right  component  replaced. 

•  Do  feature  extraction  for  Case  1  processing. 

•  If  the  replacing  component  and  the  replaced  component  are  dif¬ 
ferent,  then  find  all  instances  of  the  replaced  component  number 
in  the  active  line,  and  change  them  to  the  replacing  number. 


•  Delete  the  current  segment  and  the  segment  following  the  cur¬ 
rent  one  from  the  active  line.  Let  the  current  segment  pointer 
point  to  the  first  segment  after  the  deleted  one. 

•  Search  for  instances  of  the  terminated  component  number  in  the 
active  line.  If  there  are  no  remaining  instances,  call 
application-dependent  subroutines,  as  appropriate,  passing  the 
address  of  the  terminated  component’s  blob  descriptor. 

V  EXTRACTION  OF  FEATURES 

The  preceding  sections  outlined  a  procedure  for  isolating  connected 
components  (blobs)  from  one  another  in  an  image.  The  results  of  the 
analysis  are  a  number  of  blob  descriptor  records  that  contain  information 
about  the  blob.  Each  blob  descriptor  consists  of  several  features  or 
characteristics  of  interest. 

A.  Color 

The  simplest  feature  is  color .  The  color  of  a  blob  will  be  0 
if  the  blob  is  black  and  1  if  it  is  white.  If  the  color  of  the  background 
is  appropriately  set  up  at  initialization  time,  then  the  following  will 
ensure  that  each  blob  descriptor  subsequently  created  has  the  correct 
color:  whenever  a  new  blob  descriptor  is  created  (Case  2)  obtain  the 
color  feature  of  the  surrounding  blob  and  store  the  opposite  color  in 
the  appropriate  place  in  the  new  blob  descriptor. 

B.  Parent,  Child,  and  Sibling 

An  important  class  of  features  describes  inclusion  relationships 
among  blobs.  Because  we  have  been  careful  in  the  way  connections  are 
made  between  diagonally  adjacent  pixels,  we  can  guarantee  that  every 
component  (except  the  background)  has  one  and  only  one  surrounding 
component  blob.  We  call  any  blob’s  surrounder  its  parent,  any  holes 
(or  blobs  enclosed  by  this  one)  children,  and  other  blobs  enclosed  by 
the  same  parent  siblings. 

At  blob  creation  time  (Case  2)  a  pointer  to  the  surrounding  blob 
descriptor  can  be  stored  as  the  new  blob’s  parent  feature.  The  child 


feature  of  the  new  blob  is  set  to  0.  We  must  also  record  the  fact  that 
the  parent  blob  has  an  additional  child.  There  is  only  a  single  entry 
in  the  blob  descriptor  for  the  child  feature.  The  child  feature  of  any 
blob  points  to  the  most  recently  created  surrounded  blob;  the  sibling 
feature  of  that  most  recently  created  child  points  to  the  next  most 
recently  created  child,  and  so  forth.  The  procedure  to  follow  in  Case  2 
processing  is  as  follows:  Copy  the  child  feature  of  the  surrounding 
blob  to  the  sibling  feature  of  the  new  blob.  Set  the  child  feature  of 
the  surrounding  blob  to  point  to  the  new  blob. 

Case  1  processing  frequently  causes  two  blobs  to  be  merged  into 
one.  Topological  considerations  guarantee  that  the  two  blobs  to  be 
merged  will  have  the  same  color  feature  and  the  same  parent.  But  to 
preserve  the  integrity  of  the  parenf-child-sibling  links,  several 
actions  must  be  performed.  In  each  of  the  children  of  the  replaced  blob, 
the  parent  feature  must  be  changed  to  point  to  the  replacing  blob.  All 
the  children  of  the  replaced  blob  must  be  concatenated  onto  the  sibling 
list  of  the  children  of  the  replacing  blob.  The  pointer  to  the  replaced 
blob  itself  must  be  removed  from  the  sibling  list  that  includes  both 
the  replaced  and  the  replacing  blobs.  Finally,  if  the  common  parent  of 
the  two  merged  blobs  points  to  the  replacing  blob  as  its  child  feature, 
that  child  feature  must  be  replaced  by  the  contents  of  the  sibling 
feature  of  the  replaced  blob. 

A  greal  deal  of  bookkeeping  overhead  is  eliminated  by  recording  only 
parent  relationships.  The  children  and  siblings  can  be  regenerated 
from  the  list  of  parent  relationships.  In  the  current  implementation  of 
LINEBYLINE,  the  parent  relationship  is  not  computed  because  it  is  not 
used  by  subsequent  programs. 


C .  The  Bounding  Rectangle 

The  bounding  rectangle  of  a  blob  is  specified  by  the  minimum  and 
maximum  values  of  x  and  y  over  the  blob.  It  is  easily  extracted  as 
follows:  when  a  blob  is  created  (Case  2),  copy  the  current  value  of  y 
(which  is  constant  across  an  entire  line  of  input)  to  the  new  blob’s 
ymin  feature.  At  the  same  time  place  the  starting  column  number  (of  the 


contrasting  run-length  segment)  in  the  new  blob's  xmm  feature  and  the 
ending  column  number  in  the  blob’s  ymin  feature. 

For  each  run-length  segment  added  to  a  blob  (Case  3) ,  set  the  xmin 
feature  to  the  smaller  of  (1)  the  starting  column  number  and  (2)  the 
previous  value  of  the  xmin  feature.  Do  the  analogous  operation  for  the 
xmax  feature  and  the  ending  column  number.  The  ymax  feature  is  set  to 
the  current  line  number. 

When  two  blobs  are  merged  as  a  result  of  Case  1  processing,  the 
ymin,  xmin,  and  xmax  features  of  the  replacing  blob  must  take  the  ex¬ 
treme  values  of  the  two  blobs  to  be  merged. 

D.  Area  and  Center  of  Gravity 

The  area  of  a  blob  is  set  to  0  when  the  blob  descriptor  is  created. 

When  a  new  line  segment  of  the  blob  is  processed,  the  number  of  pixels 

in  the  segment  is  added  to  the  area  of  the  blob.  When  two  blobs  are 

merged,  their  areas  are  added  together. 

The  center  of  gravity  (CG)  is  recomputed  after  each  new  line  segment 
is  added  to  the  blob.  The  updated  CG  is  a  weighted  average  of  the  CG 
of  the  new  line  and  the  CG  of  the  previously  processed  portion  of  the 
blob.  The  weights  are  proportional  to  the  areas  of  the  new  line  segment 
and  the  previously  processed  portion,  respectively.  When  two  blobs  are 
merged,  a  weighted  average  of  the  two  CG’s  is  computed  in  a  similar 
manner  as  above. 

E.  Perimeter  Lists 

Associated  with  each  inked  blob  is  a  set  of  perimeter  points.  Be¬ 
cause  each  line  of  the  image  is  processed  sequentially,  the  perimeter 
points  are  obtained  in  a  somewhat  haphazard  manner.  During  processing, 
the  perimeter  of  a  simple  convex  blob  has  two  uncompleted  ends:  One  at 
the  starting  column  of  the  active  line  segment  (the  left  end),  and  one 
at  the  ending  column  of  the  active  line  segment  (the  right  end).  The 
perimeter  points  are  stored  as  a  singly  linked  list  called  a  perimeter 
section.  The  right  endpoint  is  at  the  beginning  of  the  list,  and  the 


B-22 


I 


left  endpoint  is  at  the  end  of  the  list.  When  there  is  more  than  one 
active  line  segment  corresponding  to  the  same  blob,  one  perimeter  section 
is  associated  with  each  active  line  segment.  The  perimeter  sections 
are  arranged  as  a  doubly  linked  ring.  This  structure  is  illustrated 
below,  using  component  2  in  the  previous  example  after  six  lines  of 
processing. 


The  structure  of  the  perimeter  section  records  and  perimeter  point 
records  in  the  illustration  above  are  shown  here: 


left _ 

right _ 

previous  section 
next  section 


Perimeter  section  record 


Perimeter  point  record 


When  each  new  line  of  pixels  is  processed,  the  lists  of  perimeter 
points  are  modified  according  to  which  case  applies.  In  all  cases 
perimeter  lists  are  compiled  for  inked  blobs  only. 

Case  3  requires  the  simplest  processing.  When  a  new  line  segment 
is  added  to  an  inked  blob,  the  left  end  of  the  associated  perimeter 
section  and  the  right  end  of  the  next  perimeter  section  are  extended. 
An  example  is  shown  below: 


Image  so  far: 


New  segment: 


Perimeter  lists 
before  new  seg¬ 
ment  is  added : 


Perimeter  lists 
after  new  seg¬ 
ment  is  added: 


bbbbbbbbbb 
b  0  1  1  1  1  1  0  0  b 
bOllllllOb 
b  0  1  1 

[0  0] 


••HE 


For  Case  1  (deletion  of  a  segment)  there  are  also  two  subcases. 
Subcase  A  occurs  if  the  deleted  segment  has  the  color  ink;  otherwise 
Subcase  B  occurs.  In  Subcase  A  two  adjacent  perimeter  sections  of  the 
terminating  blob  are  merged.  If  the  terminating  blob  has  only  one 
perimeter  section,  its  left  and  right  ends  are  joined,  and  processing 
of  this  blob  is  finished.  An  example  is  shown  below. 


Image  so  far:  bbbbbbbbbb 

b  0  0  1  1  1  1  0  0  b 
b  0  1  1  1  1  1  0  0  b 
b 

to  0] 


New  segment: 


|  In  Subcase  B  the  perimeter  sections  associated  with  the  replacing 

\  and  replaced  components  are  merged.  If  the  replacing  and  replaced 

components  belong  to  different  blobs,  the  two  blobs  will  be  merged  and 
the  two  rings  of  perimeter  sections  will  be  merged  into  one  ring.  An 
example  is  shown  below. 


Image  so  far: 


New  segment: 


bbbbbbbbbb 
bllOOOllOb 
b  1  1  0  0  1  1  0  0  b 
b 

[1111111] 


Perimeter  lists 
before  new  seg¬ 
ment  is  added: 


Perimeter  lists 
after  new  seg¬ 
ment  is  added : 


To  save  space  during  processing  of  the  image,  and  to  store  the 
completed  perimeter  lists  in  less  disk  space,  the  perimeter  lists  are 
periodically  packed.  The  variable  LINEGROUP  controls  how  often  the 
lists  are  compressed  (the  words  "pack"  and  "compress”  mean  the  same 
thing).  For  example,  when  LINEGROUP  =  25,  compression  occurs  at  lines 
25,  50,  75,  etc.  Those  points  in  the  list  that  previously  have  been 
compressed  are  skipped.  When  a  perimeter  list  has  been  completed  and 
is  ready  to  be  written  to  disk,  it  is  unpacked  and  then  repacked  if 
the  Boolean  variable  FINALPACKING  is  true.  If  FINALPACKING  is  false, 
the  perimeter  list  is  only  unpacked.  Currently,  the  value  of  FINALPACKING 
is  specified  by  the  user  at  the  beginning  of  LINEBYLXNE. 


Each  unpacked  perimeter  point  consists  of  two  16-bit  integers  with 
values  between  1  and  512.  One  integer  is  for  the  line  number  Y  and  the 


second  integer  is  for  the  column  number  X.  Compression  is  achieved  by 
converting  the  perimeter  information  from  a  series  of  (Y,X)  points  to  a 
starting  (Y,X)  point  and  a  series  of  directions  from  one  point  to  the 
next.  The  starting  point  has  the  same  format  as  an  unpacked  perimeter 
point,  except  that  512  is  added  to  the  column  number.  Therefore,  a 
column  number  greater  than  512  identifies  that  point  as  the  beginning  of 
a  series  of  compressed  points.  Because  the  perimeter  is  8-connected,  there 
are  eight  possible  directions  indicating  the  positions  of  subsequent 
perimeter  points.  They  are  numbered  from  0  to  7  and  are  assigned  as 
follows : 


Change  in  X 
-10  1 

-13  2  1 

Change  in  Y  0  4  0 

15  6  7 

For  example,  if  the  perimeter  list  of  the  last  example  above  (deletion 
of  a  segment,  Subcase  B)  were  to  be  compressed  starting  with  the  point 
(3,1),  the  direction  to  the  next  point  (2,1)  would  be  2;  the  direction 
from  (2,1)  to  the  next  point  (1,1)  would  be  again  2;  and  the  direction 

from  (1,1)  to  the  next  point  (1,2)  would  be  0. 

Each  of  the  eight  directions  can  be  represented  as  a  3-bit  code 
(000  *  0,  001  *  1,  010  *  2,  ...,  Ill  *  7).  Five  direction  codes  fill 
15  bits  of  a  16-bit  word.  The  16th  bit  (the  sign  bit)  is  set  equal  to  1 
to  Indicate  that  the  word  contains  compressed  perimeter  points.  Each 
set  of  five  direction  codes  is  stored  alternately  in  the  line  and  column 
locations  of  the  perimeter  point  record.  If  we  continue  the  example, 
the  first  five  directions  would  be  2,  2,  0,  6,  and  7.  The  would  be 
coded  and  placed  in  a  16-bit  word  as: 

1  111  110  000  010  010 

16  15  14  13  12  11  10  9  8  7  6  5  4  3  2  1 
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This  word  is  interpreted  as  the  two’s  complement  integer  -1756  octal, 
or  -1006  decimal.  The  second  five  directions  are  0,  1,  1,  0,  and  5. 
Their  bit  representation  is: 

1  101  000  001  001  000 

16  15  14  13  12  11  10  9  8  7  6  5  4  3  2  1 

This  is  interpreted  as  -27670  octal,  or  -12216  decimal.  Hence,  the 
first  11  perimeter  points  are  compressed  into  the  space  of  two  perimeter 
point  records  as  follows:  (3,513),  (-1006,-12216). 

Compression  of  groups  of  less  than  10  points  is  not  attempted. 
Therefore,  during  processing,  the  list  of  compressed  points  is  inter¬ 
spersed  with  uncompressed  points.  In  addition,  there  are  usually  a 
few  uncompressed  points  at  the  end  of  the  list.  In  our  example,  the 
12th  point,  (3,8)  remains  uncompressed.  The  completely  compressed 
list  is  therefore: 


The  unpacking  process  is  just  the  reverse  of  the  packing  process. 
The  column  number  of  each  perimeter  point  is  examined.  If  it  is  greater 
than  512,  the  start  of  a  series  of  compressed  points  has  been  detected. 
Subtracting  512  from  the  column  number  results  in  the  (Y,X)  coordinate 
of  the  starting  point.  The  sign  of  the  column  number  of  the  next 
perimeter  point  is  examined.  If  it  is  less  than  or  equal  to  0,  then 
compression  of  10  points  has  been  detected.  The  locations  of  these 
points  are  recomputed  from  the  3-bit  direction  codes  that  are  extracted 
from  the  two  16-bit  words.  This  process  continues  until  a  perimeter 
point  record  with  a  column  number  greater  than  0  is  encountered. 


VI  HIERARCHY  OF  PROCEDURES 


This  chart  shows  the  hierarchy  of  procedures  within  LINEBYLINE. 
All  of  the  procedures  that  are  called  by  each  procedure  are  indented 
and  listed  below. 


LINEBYLINE 

PACKPERIMS 

PKPERIMS 

COMPRESSED 

UNPACKPERIMS 

UNPACK5 

DIRPOINT 

DELETEPERIMS 

NEWPERIM 

DIRECTION 

GETLINE 

ENDOFBLOCK 

SETLINE 

CREATERUNLENGTHS 

ADDRUNLENGTHS 

PROCESSLINE 

DELETESEGMENT 

ADDLEFTPERIMPOINT 

NEWPERIM 

ADDRIGHTPERIMPOINT 

NEWPERIM 

RECORDBLOB 

WRITEPERIMS 

WRITEPOINTS 

PKPERIMS 

(see  above) 
COUNTANDWRITE 
COMPRESSED 
UNPACKPERIMS 

(see  above) 
WRITEPOINTS 

(see  above) 
DELETEPERIMS 

WRITEBLOB 

DELETEBLOB 

INSERTSEGMENT 

NEWBLOB 

ADDRIGHTPERIMPOINT 
(see  above) 

CASE3PROCESSING 

UPDATEPERIMS 

ADDLEFTPERIMPOINT 
(see  above) 
ADDRIGHTPERIMPOINT 
(see  above) 
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VII  INPUT  IMAGE  FORMAT 


The  binary  image  that  is  read  by  LINEBYLINE  is  512  *  512  pixels. 
There  are  8  pixels  per  byte  (each  pixel  is  one  bit).  Therefore,  one 
line  of  data  occupies  64  bytes,  and  the  entire  image  occupies  512  *  64  = 
32K  bytes  on  disk.  The  data  is  read  one  byte  (one  character  in  PASCAL) 
at  a  time. 

However,  there  is  one  small  difference  due  to  the  PDP  11/40  version 
of  PASCAL.  The  bytes  are  stored  on  disk  in  blocks  of  512  bytes.  The 
PASCAL  I/O  processor  interprets  a  block  boundary  as  an  extra  byte. 

To  compensate  for  this  inaccuracy,  the  LINEBYLINE  software  skips  an 
"imaginary"  byte  once  every  512  times  (when  the  Boolean  function 
ENDOFBLOCK  is  true).  On  most  computer  systems,  this  would  not  be  neces- 
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ANNOTATIONS  FOR  THE  SPSS  RESULT  PACKAGE 
(Submitted  Separately) 


The  printout  of  an  SPSS  analysis  contains  several  parts*  Following 
is  a  description  of  each  part  with  the  page  number  of  the  listing  on 
which  that  part  of  the  printout  begins.  Preceding  the  actual  output, 
some  initial  information  is  listed,  including  the  control  file  used 
to  run  SPSS  and  the  data  that  are  to  be  read  by  SPSS.  These  data  are 
ordered  as  a  sequence  of  trials.  The  numbers  represent:  Subject  code 
(1-30),  session  number,  trial  number,  and  29  feature  values. 

•  Page  1:  SPSS  description;  format  of  the  data  that  are  to  be 
read  in. 

•  Page  3:  A  partial  list  of  the  data  that  were  read  in.  Each 
trial  is  listed,  but  only  five  selected  feature  values  are 
listed. 

•  Page  9:  SPSS  options  used. 

•  Page  10;  Number  of  trials  per  subject  (group). 

•  Page  10 :  Feature-value  means  for  each  subject. 

•  Page  12 :  Feature-value  standard  deviations  for  each  subject. 

•  Page  15 :  Wilkfs  lambda  and  univariate  F-ratios  for  each  feature. 

Each  of  these  numbers  is  related  to  the  discriminating  power  of 
the  specific  feature.  For  Wilk's  lambda,  the  smaller  the 
number,  the  better  is  the  discriminating  power  of  that  feature. 

The  opposite  is  true  for  F-ratios:  The  larger  the  F-ratio,  the 
greater  is  the  discriminating  power  of  that  feature. 

•  Page  17:  Within-groups  correlation  matrix.  This  is  similar  to 
the  covariance  matrix,  except  that  each  feature  mean  is  first 
normalized  to  mean  *  0  and  standard  deviation  =  1  before  calcu¬ 
lation  of  the  covariance.  A  number  in  this  matrix  that  approaches 
1  indicates  a  high  degree  of  correlation  between  those  two 
features.. 

/ 

0  Page  19:  Prior  probability  of  finding  any  particular  subject 
within  the  population.  A  priori,  this  probability  was  set  to 
be  equal  for  each  subject. 
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•  Page  1^:  Twenty-nine  discriminant  taiu  tio;u  w<r.i determined. 

The  data  in  this  chart  demonstrate  the  declining  discriminating 
power  of  each  additional  discriminant  function. 

*  Page  20:  The  standard  discriminant-function  coefficients  are 
listed.  The  product  of  all  normalized  feature  values  and  the 
appropriate  discriminant-function  coefficient  will  be  the  dis¬ 
criminant  score  used  to  classify  each  subject. 

•  Page  23:  Here  start  listings  of  both  the  unstandardized 
discriminant-function  coefficients  and  the  location  in  29- 
dimensional  space  of  the  centroids  of  each  group.  The  centroid 
is  the  mean  discriminant-function  score  for  each  subject. 

#  Page  30:  Classification  results:  Actual  group  is  the  subject 
code  number;  "highest  group  D**2"  is  the  classification  result. 

In  all  cases  the  group  was  classified  correctly.  "Probability 
P(G/x)"  is  the  probability  that  the  sample  belongs  to  Group  G 
given  that  the  discriminant  measure  x  was  observed.  "Probability 
P(x/G)"  is  the  probability  that  the  discriminant  measure  x  would 
be  observed  given  that  the  sample  is  from  Group  G.  The  actual 
discriminant  scores  (upon  which  the  classification  is  based)  are 
listed  on  the  far  right. 

#  Page  78:  Territorial  map  of  first  two  (most  discriminating) 
discriminant  functions.  This  is  a  2D  representation  oi:  rh'>  data. 
The  centroids  (*)  are  separated  by  numbers  or  letters  (correspond- 
ing  to  that  particular  group,  numbered  0-9  or  A-Z)  that  define 
the  equal  probability  lines  between  groups. 
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MISSION 

of 

Ram  Air  Development  Center 

RA DC  plans  and  executes  research,  development,  tut  and 
selected  acquisition  programs  In  support  oi  Command,  Control 
Communications  and  Intelligence  (C^I)  activities.  Technical 
and  engineering  support  within  areas  oi  technical  competence 
Is  provided  to  BSD  Program  OHices  (PCs)  and  other  BSD 
elements.  The  principal  technical  mission  areas  are 
communications ,  electromagnetic  guidance  and  control,  sur¬ 
veillance  oi  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  iniormation  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


